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Abstract 

The  combined  capacity  expansion,  production- 
inventory  problem  Is  modeled  as  a  linear,  integer 
program  and  then  extended  to  also  allow  capacity  de- 
crease.  The  models  assume  constant  returns  to  scale 
in  the  production  function  of  a  firm  which  must  meet 
deterministic  demands  for  a  single  product  over  N 
periods,  with  no  backordering,  at  minimum  cost. 

A  linear  transformation  is  used  to  obtain 
equivalent  forms  of  the  models  which  are  then  decompos- 
ed into  fixed  cost  and  variable  cost  parts.   Enumer- 
ation is  done  on  the  fixed  cost  variables  and  transpor- 
tation sub-problems  are  solved  with  the  remaining  var- 
iables.  Special  demand  and  cost  structures  are  inves- 
tigated, one  of  which  leads  to  formulation  of  the  en- 
tire problem  as  a  network. 

Extensions  are  made  to  include  piecewise- 
linear,  convex  production  costs  and  lumpiness  in  ca- 
pacity acquisition.   A  reinterpretation  of  the  models 
to  handle  the  employment  smoothing  problem  is  also  pre- 
sented.  Finally,  computational  results  of  a  Fortran  IV 
coding  of  the  basic  algorithms  are   presented. 
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Chapter  I 
I. 1  Introduction 

The  quantitative  Rpproach  to  the  solution 
of  problems  In  the  realm  of  managerial  economics  can 
be  traced  back  to  the  efforts  of  the  early  Industrial 
engineers  In  the  United  States  at  the  end  of  the  nine- 
teenth century.   These  men,  most  notably  Frederick  W. 
Taylor,  pioneered  the  Industrial  engineering  concepts 
such  as  Job  description,  time-motion  studies,  etc., 
which  had  such  a  great  Influence  on  the  development 
of  the  assembly  line  In  Industry  during  the  first 
half  of  the  twentieth  century. 

With  World  War  II,  however,  an  entirely  new 
set  of  tools  was  refined  and  applied  by  the  early  prac- 
titioners In  the  Infant  field  of  Operations  Research. 
Though  most  of  the  techniques  had  their  origin  In 
military  applications,  they  gradually  found  their  way 
Into  the  Industrial  sphere  In  the  peace-time  economic 
expansion  resulting  from  the  war.   In  the  Intervening 
decades,  the  complexity  of  the  practical  problems 
solved  has  moved  apace  with  the  rapid  advances  In  theory 
typical  of  a  new  field.   The  logical  result  is  that  an 
increasing  amount  of  effort  has  been  put  into  inte- 
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2. 
gratlon  of  past  techniques  using  newer  theory  in  an 
attempt  to  provide  more  realistic  approaches  to  problems 
with  broader  horizons. 

The  present  thesis  makes  such  an  attempt  by 
offering  models  of  the  combined  production-inventory 
and  capacity  expansion  problem  and  algorithms  which 
give  all  integer  solutions  to  these  models.   Dis- 
cussion of  the  importance  of  the  combined  problem  and 
the  integral  approach  are  given  later  in  this  chapter. 
This  is  the  first  integer  approach  to  the  combined 
problem  known  to  exist. 

The  basic  technique  used  will  be  branch  and 
bound  with  network  sub-problems,   A  unique  feature 
of  the  algorithm  is  the  capability  of  starting  it 
at  any  point  in  the  enumeration.   This  is  especially 
meaningful  when  linear  programming  is  uBed  to  give  a 
universal,  and  often  non-integer,  optimal  solution. 
The  enumeration  can  be  started  from  a  suitable  point 
with  respect  to  the  linear  programming  optimum  to 
determine  the  existence  of  any  optimal,  rounded-off 
solutions  rather  quickly. 

The  famous  economic  problem  of  "lumpiness" 
in  capacity  acquisition  as  well  as  multi-shift,  and 
overtime,  operation  of  the  production  facility  will 
be  treated,   A  discussion  of  the  relationship  between 
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the  traditional  capacity  expansion  and  employment 

smoothing  literatures  and  re-interpretation  of  the  re- 
sults in  terms  of  employment  smoothing  will  also  be 
given, 

1,2   Problem  Description 

A  useful  paradigm  is  a  company  with  a  single 
producing  facility  and  a  homogeneous  product  faced 
with  a  schedule  of  deterministic  demands  over  a  set  of 
discrete  periods  to  some  pre-set  horizon.   These  de- 
mands must  be  met  on_  time  from  production  or  existing 
inventory  at  minimum  possible  cost. 

The  capacity  of  the  production  facility  is 
defined  as  its  output  in  units  per  time  period  when 
equipment  is  being  run  at  its  maximum  efficient  rate 
with  a  single  shift  not  working  overtime.   The  models 
will  be  extended  to  allow  multiple  shifts  and  over- 
time work,  but  note  that  these  additions  are  not  ger- 
mane to  the  basic  definition  of  capacity. 

The  company's  basic  problem,  assuming  that 
cumulative  demand  never  exceeds  cumulative  capacity 
(where  cumulative  capacity  now  includes  all  regular 
and  overtime  production),  is  to  find  a  production- 
inventory  schedule  to  meet  demand  at  minimum 
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cost.^  References  to  this  problem  are  ubiquitious  In 
both  the  operations  research  literature  and  the  ec- 
onomics literature,  and  will  not  be  discussed  in  this 
work. 

If,  however,  cumulative  demand  does  exceed 
cumulative  capacity  within  the  planning  horizon,  then 
the  company  must  make  a  capacity  expansion  decision. 
That  is,  it  must  decide  at  what  time  Cor  times)  and  by 
what  amount  Cor  amounts)  to  invest  in  capital  equipment 
so  that  production  can  satisfy  demand.   The  intuitive 
difficulty  of  this  problem  is  the  fact  that  a  change 
in  the  capacity  in  one  period  affects  the  production- 
inventory  schedule  for  all  future  periods,  thus  lending 
combinatorial  aspects  to  the  capacity  expansion  problem. 
Further  difficulty  arises  when  capacity  must  be  added 
in  indivisible  units,  each  capable  of  producing  large 
numbers  of  units  of  output  as  is  usually  the  case  in 
industry. 

Capital  equipment  itself  has  many  aspects. 
Deterioration  might  be  a  significant  factor  with  the 


-'•The  existence  of  capacity  in  excess  of  demand 
does  not  preclude  a  decision  to  expand  capacity  since 
there  could  exist  circumstances  where  such  a  move 
would  lead  to  lower  costs. 
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type  of  capital  Installed,  in  which  case  the  capacity 

expansion  problem  is  further  complicated  by  the  need 
for  equipment  replacement.   On  the  other  hand,  fixed 
capital  might  have  a  relatively  long  life,  making  per- 
iodic  maintenance  the  only  significant  factor  within 
the  time-span  of  interest.   It  is  also  possible  that  the 
fixed  capital  equipment  can  be  salvaged  or  scrapped  if 
warranted  by  changes  in  the  demand  situation.   The  first 
case  can  only  be  approximated  by  the  models  which 
follow  and  then  only  by  allowing  maintenance  costs  to 
reflect  the  increased  investment  over  time  needed  to 
maintain  equipment  at  its  original  capacity.   In  add- 
ition to  not  obtaining  equipment  replacement  decisions 
from  the  models,  the  aj^proximation  is  faulty  in  that 
the  increased  maintenance  costs  are  applied  uniformly 
to  equipment  in  each  period,  regardless  of  age.   The 
latter  two  cases  are  modeled  exactly,  however,  the  sec- 
ond by  the  explicit  structure  of  the  maintenance  costs 
gind  the  last  implicitly  by  the  relevsmt  mathematical 
programming  model, 

I. 3  Economic  Background 

With  respect  to  the  economic  geography,  the 
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problem  lies  on  the  boundary  of  capital  theory  and  the 
classical  theory  of  production.   A  good  discussion  of 
the  problems  arising  from  the  merger  of  these  two  dis- 
ciplines is  contained  in  the  unpublished  Ph.  D, thesis 
of  Carroll  (Carroll  68). 

Apropos  of  capital  theory,  a  decision  must 
be  made  as  to  the  proper  points  in  time,  and  the  proper 
amounts,  for  investment  in  capacity;  however,  there  is 
no  requirement  present  for  alternative  investments. 
In  general,  if  expansion  of  capacity  in  a  number  of 
different  areas  were  being  considered  in  each  period 
with  a  limited  budget,  then  the  extensive  work  in  cap- 
ital budgeting  would  be  relevent.^  Further,  capital 
budgeting  problems  typically  abstract  from  the  micro 
problem  of  optimizing  the  form  of  each  investment  pro- 
ject, being  concerned  mainly  with  the  macro  allocation 
of  funds  to  projects.   The  micro-optimization  is  an  im- 
portant part  of  the  work  that  follows.   Extensions  of 
the  present  problem  to  the  capital  budgeting  sphere 
will  be  discussed  in  the  final  chapter. 


^A  classical  application  of  mathematical  program- 
ming to  the  capital  budgeting  problem  is  given  by 
Weingartner  (Weingartner  63). 
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In  the  theory  of  production  as  dlacussecl  by 
Carroll,  consideration  of  production.  Inventory,  and 
capacity  expansion  places  the  problem  In  between  the 
short  term,  where  capacity  Is  fixed  and  decisions  are 
only  made  on  levels  of  production,  and  the  long-run  de- 
cision of  optimal  overall  capacity  level,  which  ab- 
stracts from  all  short-run  decisions.   In  effect,  the 
mechanism  of  inventorying  production  functions  as  the 
mediator  between  the  short-run  minimum  cost  production 
problems  and  the  long-run  minimum  cost  capacity  expan- 
sion problems,  serving  to  "smooth"  the  extremes  in  the 
optimal  solutions  to  the  individual  problems. 

A  more  satisfying  concept  of  the  spectrum  of 
decision  time-spans  is  given  by  Maxwell  (Maxwell  65). 
He  considers  the  very  short-run  as  the  time-span  in 
which  all  factors  of  production  are  fixed  in  both  ac- 
quisition and  use.   Successive  decision  time-spams  are 
encountered  as  factors  become  variable  in  use  and  ac- 
quisition, giving  a  succession  of  short  runs  of  in- 
creasing length.   The  long-run  is  finally  achieved  when 
all  factors  become  variable  in  both  use  and  acquisition. 
Against  this  spectrum,  the  problem  at  hajid  ranges  from 
a  succession  of  relatively  long  short-runs  when  only 
capital  stock  is  fixed  in  acquisition  and  all  factors 
are  variable  in  use,  to  a  long-run  over  which  the 
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capital  stock  is  variable  in  acquisition  at  a  succession 
of  discrete  points.   Correspondence  with  classical 
theory  is  not  exact  though,  since  capital  is  durable 
and  not  perfectly  flexible  as  required  in  the  classical 
long-run.   The  necessity  for  providing  inter-period  rela- 
tionships  through  the  mechanism  of  inventorying  in  the 
case  of  inflexible  capital  was  pointed  out  by  Carroll: 


...So  long  as  the  theory  of  production  ab- 
stracted from  the  particular  problems  of  fixed 
In-puts,  it  was  essentially  a  mono-period 
theory  in  which  inter-period  relationships 
could  largely  be  neglected.   Even  the  long- 
run  analyses  tended  to  treat  the  question  of 
determining  the  optimal  capital  stock  as  a 
once-and-for-all  decision.   Since  capital  was 
generally  assumed  to  be  completely  flexible 
in  the  long-run,  questions  of  the  Impact  of 
one  period's  capital  accumulation  decision 
upon  that  of  another  were  not  raised.   But 
once  capital  came  to  be  viewed  as  a  durable 
good,  which  was  not  perfectly  flexible, 
inter-period  (i.e.,  between  long-run  periods) 
relationships  were  a  matter  of  Interest. 
(1968,  p.  5) 


In  the  following  work,  a  salvage  value  of  capital 
goods  is  Introduced  and  constrained  to  be  less  than  the 
acquisition  cost,  thus  ensuring  at  least  a  mild  flex- 
ibility in  capital. 

Maxwell  also  presents  a  variation  on  the  con- 
ventional theory  of  the  firm  which  is  relevant  here. 
The  contrast  is  in  the  derivation  of  the  production 
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function  which,  in  the  classical  theory,  is  obtained 
from  the  combination  of  varying  inputs  of  consumptive 
factors  of  production  with  the  fixed  factor  to  obtain 
a  graph  of  maximum  outputs  over  the  range  of  the  con- 
sumptive factors.   The  underlying  assumption  is,  of 
course,  that  the  fixed  factor  is  indivisible  in  use  and 
must  be  used  to  the  available  limit.   If  instead,  the 
fixed  factor  is  viewed  as  indivisible  in  acquisition, 
but  divisible  in  use,  a  more  useful  approach  to  the 
firm's  production  function  can  be  derived. 

Assume  constant  returns  to  scale  in  the  pro- 
duction function.   Then,  for  each  level  of  input  of 
consumptive  factors,  the  fixed  factor  is  used  up  to  its 
zero  level  of  marginal  productivity  giving  a  linear  ex- 
pansion path  for  the  firm  up  to  the  capacity  of  the 
fixed  factor.   Thus,  under  constant  returns,  the  scale 
of  use  of  the  fixed  and  consumptive  factors  of  produc- 
tion is  varied  linearly  up  to  the  capacity  of  the  fix- 
ed factor,  while  the  proportions  of  both  remain  cons- 
tant.  Beyond  the  capacity  of  the  fixed  factor  the  pro- 
portion of  the  consumptive  factor  is  increased,  giving 
returns  to  proportion  in  production  as  opposed  to  the 
returns  to  scale  obtained  in  the  former  range.   The 
importance  of  this  concept  to  the  theory  of  the  firm  is 
pointed  out  by  Maxwell: 


•r  •!         « 
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Many  empirical  cost  studies  indicate  that  the 
marginal  costs  tend  to  be  constant  very  nearly 
up  to  capacity,  suggesting  that  this  range  in 
which  less  than  all  of  the  stock  of  fixed- 
factor  services  is  being  utilized  may  indeed 
be  a  significant  portion  of  the  total  range 
of  outputs  that  the  firm  is  capable  of  pro- 
ducing in  the  short-run,   (1965,  p.  lo6) 


For  the  problem  of  interest,  the  range  of  the 
production  function  reflecting  returns  to  proportion  is 
of  no  interest  since  one  of  the  decision  variables  for 
each  period  is  the  optimal  level  of  fixed  capital  for 
that  period.   Thus,  concern  can  be  focused  solely  on  the 
constant  returns  to  scale  portion  of  the  production 
function  assumed  present  up  to  the  capacity  of  the  fix- 
ed factor. 

In  the  problem  which  follows,  the  bundle  of 
consumptive  factors  of  production  9^  which  is  needed 
in  period  i  for  one  unit  of  output  is  assumed  known. 
Further,  it  is  assumed  that  the  fixed  factor  of  pro- 
duction is  used  at  its  zero  marginal  productivity  level 
for  this  bundle  d^   or  any  multiple  of  it.   The  capacity 
W.  of  the  fixed  factor  in  any  period  can  then  also  be 
expressed  in  terms  of  the  maximum  number  of  bundles  3. 
which  can  be  used  efficiently  with  the  fixed  factor  in 
the  period  (in  this  case  the  period  includes  only  reg- 
ular time  work  by  one  shift).   The  level  of  production 
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in  period  1  which  is  x^  must  thus  satisfy 

X.  -  W^     1  =  1,.. . ,N 

These  assumptions,  when  taken  in  the  context 
of  fixed  capital  investment,  suppress  the  problem  of 
labor  force.   That  is,  although  the  labor  force  is  a 
non-consumptive  factor  of  production,  the  necessary 
assumption  here  is  that  labor  force  is  always  adjusted 
to  the  exact  size  required  to  run  the  capital  equipment 
at  the  predetermined  optimal  level  x^^  in  period  i.   An 
interpretation  of  the  problem  of  capacity  in  terms  of 
labor  force  rather  than  capital  equipment  will  be 
given  subsequently. 

The  objective  is  to  minimize  the  present 
value  of  the  costs  incurred  over  N  periods  to  meet  a 
predetermined  demand  schedule.   It  should  be  noticed 
that  the  problem  is  abstracted  from  pricing  consider- 
ations and  from  determination  of  optimum  output 
schedules;  in  fact,  concern  is  with  the  optimum  method 
of  meeting  any  demand  schedule.   To  achieve  the  pre- 
sent value  criterion,  it  is  only  necessary  that  the 
individual  costs  in  the  objective  function  of  the  model 
be  discounted  by  the  appropriate  cost  of  capital  p^^  in 
period  i.   Thus  discounting  does  not  have  to  be  modeled 
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but  can  be  Introduced  explicitly  in  the  relevant  model 
coefficients.   The  subject  of  the  proper  choice  of  the 
discount  factor  p  is  not  germane  to  this  thesis  and 
will  not  be  discussed. 

The  major  stumbling  block  to  the  merger  of 
the  capital  and  production  theories  has  historically 
been  the  problem  of  fixed  capital  as  a  non-consumptive 
factor  of  production.   This  question  is  also  lucidly 
discussed  by  Carroll.   To  summarize,  some  method  must 
be  found  to  measure  the  value  of  a  factor  of  production, 
i.e*  fixed  capital,  which  is  not  consumed  in  the  pro- 
ductive process,  but  which  must  be  present  for  produc- 
tion to  take  place.   In  this  thesis,  the  value  measure 
will  be  the  reduction  in  cost  obtained  in  the  produc- 
tion of  a  given  demand  schedule,  a  measure  which  is  im- 
plicit in  mathematical  programming  optimization. 

Finally  in  the  economic  context,  mention 
should  be  made  of  the  historic  bugbear  of  marginal 
analysis,  indivisibility  in  acquisition  of  capital 
stock.   This  discreteness  in  acquisition  possibilities 
destroys  the  differentiability  of  most  useful  functions 
at  best,  and  their  continuity  itself  at  worst.   As  will 
be  described  later,  a  simple  adjustment  in  the  integer 
programming  formulation  will  handle  this  problem 
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fairly  easily  and,  in  some  cases,  even  improve  the  al- 
gorithms. 

I. h     Problem  Formulation 

The  interval  of  time  under  consideration  is 
divided  into  n  periods  of  equal  lengths  with  the  first 
period  designated  as  1  and   the  n^h^  or  horizon,  period 
designated  N,   It  is  assumed  that  demand  in  period  i, 
dj^,  is  knowi  deterministically  for  i  =  1,..,,n.   Then 
using  the  following  variables: 

x^-  production  in  period  i,  i  =  1,.,.,N 
v.-  inventory  from  period  i  to  period  i  +  1 
w^-   capacity  added  in  period  i,  i  =  1,..,,N 
the  following  coefficients: 

p^-   cost  of  a  unit  bundle  3^of  consumptive 
factors  of  production  in  period  i, 
i  »  1 , . . .  ,N 
h^-  unit  inventory  cost  from  period  i  to  1  +  1 
Cj^-   construction  cost  of  a  unit  of  capacity 

in  period  i 
m^-  maintenance  cost  of  a  unit  of  capacity  in 
period  i 

where  coefficients  are  assumed  integer  and 
p^  •  0,  hj_  -  0,  Ci  -  0,  m^  -  0   i  =  1,...,N 
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and   the    following  explicit    mtep;er  variables: 

Wq-  Initial   capacity  (w     -  0) 

Vq-  Initial   Inventory  (v     -  0) 

If  all   functions   are   assumed  linear,   the   capacity  ex- 
pansion-production-inventory problem   (CPI)    can  be 
written  as   the   following  Integer  programming  problem: 


N   ^  1 

PI:        min   Z^  »      2    [piXi+hiVjL+CiWj^+mi*  (    Z     wj+w^pj 


s,  t. 

(Pl.l)  X4    -      E  w.   ^  w^  1   -    1,...,N 

1        J=l  JO  »        » 

x^   -   vi   -    d^  -    Vq 
CPI. 2)        Xi  +    v^_2.  -    v^   »   di  1   »    2,...,N-1 

X,    V,   w  -   0        and  INTEGER 


where   the  right-hand  side   of   the   constraints   is   assumed 
integer.      Constraints    (Pl.l)   restrict  production   to  be 
less   than   or  equal  to   capacity   in   any  period,   and  con- 
straints   (Pl,2)    require   demand  to   be  met  exactly   out 
of  the   production  and   inventory   in   each  period.      All 
costs   are   non-stationary   and  are   explicitly   discounted 
to  the   first  jperiod  by   the  appropriate  discount   rate  pj^, 
which  is   external   to  the  model  as   discussed  in  section 
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1.3. 

Initial  Inventory,  Vq,  and  Initial  capacity, 
Wq,  are  set  explicitly  to  the  values  relevant  to  the 
problem  being  solved.   Initial  Inventory  and  Initial 
capacity  are  rather  straightforward  In  the  problem  con- 
text and  no  difficulty  In  obtaining  them  should  occur 
In  the  real  world.   Pinal  Inventory,  however.  Is  more 
of  a  subjective  variable.   If  Vjq  Is  allowed  to  be  an 
Implicit  variable,  then  Vj,  =»  0,  since  If  v^,  >  0  then 
the  level  of  production  In  periods  N,  N-1,  ...  could 
be  reduced  until  v„  ■  0,  giving  a  lower  cost  while 
still  minimizing  cost  (recall  that  cost  coefficients 
are  non-negative).   On  the  other  hand,  v^j  can  be  use- 
ful as  an  explicit  variable  as  shown  by  an  example. 
Assume,  for  N  *  3,  that  the  demand  schedule  Is 
d*^  «  (1,7,20).^  In  this  case,  v^   could  subjectively 
be  set  to  some  value  greater  than  zero  as  a  surrogate 
for  an  expectation  of  continued  Increase  In  demand  be- 
yond the  horizon. 

The  objective  function  In  PI  can  be  re- 


3ln  the  remainder  of  this  thesis,  an  untransposed 
vector  Is  assumed  to  be  a  column  vector. 
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written  as 


mln  Zt  -   E  [pi  X4+hiVi  +  (ci+  L   in.  )w.  ]  +  w^*  i   m. 
J-    1-1   -L  1   J.  -L    1  j.^  J   1      °  1-1  i 


It  can  be  clearly  seen  from  this  aggregated  function, 
that  the  total  cost  of  an  Incremental  unit  of  capacity 
wi  Is  equal  to  Its  acquisition  cost  In  period  1  plus 
the  discounted  sum  of  the  maintenance  costs  for  one 
unit  of  capacity  up  to  the  horizon.   Maintenance  costs 
themselves  are  period  dependent  (  I.e.,  Independent  of 
the  age  of  the  Installed  capacity). 

The  CPI  model  can  be  extended  easily  to  a 
capacity  Increase-decrease,  production-Inventory  model 
(CIDPI)  as  shown  below.   Let: 

ri-  capacity  scrapped  In  period  1 
f  j^-  salvage  value  per  unit  of  capacity 
scrapped  In  period  1 
all  other  variables  and  coefficients  are  as  before  and 
assume 

f^  <  c^     1  »  1,...,N 

fj  +  ^I  mj,  <  ci  +  Jj^mk    V  J  >  1,  1  -  1,...,N-1 

Then  the  CIDPI  cajn  be  written  as  the  following  Integer 
programming  problem: 


•.  ..  -1 
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P2:    mln  Z^   «  .^  ,  LPi^l+hiVj^+CiWi-f^rj^+mi' (  t  (wj-Pj  )+Wq)] 


S.t, 

1 

(P2.1)    x^  -  ^E^(wj-rj)  ^  Wq     1=1,...,N 

XI  -  vi  -  di  -  Vq 
(P2.2)    Xi  +  Vi_i  -  Vi  =  d^     i  =  2,...,N-1 
^N  +  vj^.i  =  dj^  +  vj^ 

i         > 

(P2.3)    Wq  +   I  (wj-rj)  -  0     1  »  1,...,N 

X,  V,  w,  ?  -  0    and  INTEGER 

where  the  right-hand  sides  of  the  constraints  are  in- 
teger as  before,  amd  constraints  (P2,l)  and  (P2.2)  have 
the  same  interpretation  as  previously.   Constraints 
(P2.3)  ensure  that  current  capacity  is  always  non- 
negative. 

It  should  be  noted  here  that   P2   contains 
the  flexibility  of  allowing  f^  <  0  as  might  be  the  case 
if,  instead  of  a  salvage  value  for  equipment,  there 
were  a  positive  cost  of  scrapping  capacity  in  some 
period  (in  which  case  scrapping  might  be  a  means  of 
avoiding  unnecessary  maintenance  cost  on  unneeded 
equipment).   For  f.  -  0,   P2   is  bounded  because  of  the 
initial  assumptions  made  on  c^^  and  f^,  and  for  fj^  <  0, 
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P2      Is   bounded  because    of   the    form  of  the    objective 
function. 

As   with      PI    ,    the   objective    function   of      P2 
can   be   aggregated  as    follows: 


mln   Z2   =      5,CPlXi+h^Vj^+(Ci+   Z   m.  )wj^   -    (fi+   I   inj)rj] 


N  N  N 

^   E    m.)w.    -    (fi+    ; 
J  =  l    J      ^  ^    J  = 

N 

+  w^*    Z   m-i 
o   1=1   1 

PI  becomes  a  special  case  of  P2   If  the 
scrapping  cost  Is  set  to  Infinity  (I.e.,  T^=   -ooVi)  in 
P2  .   Then  r  =  0  always  and  constraints  (P2.3)  become 
redundant,  giving  back  the  form  of  PI  .   Two  distinct 
models  will  be  retained,  however,  for  both  economic  and 
theoretical  reasons. 

Finally,  in  the   CIDPI  model,  the  value  of 
fixed  capital  is  central  to  the  minimum  cost  solution. 
During  periods  1  to  N,  this  value  is  reflected  in  the 
salvage  coefficient  f.   To  be  complete,  though,  the 
model  must  include  the  value  of  fixed  capital  after 
period  N,  and  here  two  cases  are  relevant;  either  the 
firm  expects  to  go  out  of  business  after  period  N  or  it 
expects  to  continue  production,  but  its  planning  horizon 
is  limited  by  dependable  demand  projections.   In  both 
cases,  a  dummy  period  N  +  1  can  be  added  to  the  model 
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with  dj^4.jL  =  0,   Then,  in  the  former  case,  fN+l  Is  Just 
the  relevant  salvage  value  of  the  fixed  capital,  while 
in  the  latter  case,  fj^^.-,  can  be  set  to  equal  some  dis- 
counted estimate  of  the  future  value  of  production 
possible  with  the  installed  capacity. 

As  for  the  CPI  model,  nothing  need  be  done 
since  the  salvage  value  of  capital  has  no  bearing  on 
the  problem.   That  is,  capital  only  has  value  in  the 
problem  context  to  the  extent  that  it  can  meet  the  de- 
mands. 

I. 5  Related  Literature 

A  large  literature  exists  in  the  area  of  pro- 
duction planning  and  there  is  at  least  a  modest  one  on 
capacity  expansion.   But  the  combined  problem,  though 
economically  significant,  has  been  given  relatively 
little  attention,  and  no  attention  at  all  has  been 
given  to  the  necessity  of  providing  good  integer  solu- 
tions  to  such  a  problem,   A  related  literature  closely 
allied  with  the  combined  problem  does  exist,  however, 
in  the  area  of  employment,  or  work- force,  smoothing. 
Each  of  these  areas  will  be  discussed  subsequently. 

Mention  must  first  be  made  of  the  relation- 
ship among  the  existing  capacity  expanaion  literature 
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economics,  and  reality.   Any  demand-related  problem 
naturally  depends  on  a  reasonable  projection  of  demand, 
either  stochastic  or  deterministic,  within  a  believable 
horizon.   Since  forecasting  is,  at  best,  a  risky  busi- 
ness, the  veracity  of  any  demand  projection  must  have 
an  inverse  relationship  to  the  length  of  the  horizon 
over  which  the  projection  is  made.   A  reasonable  figure 
for  most  companies  would  probably  not  exceed  five  years. 
Within  such  a  time  span,  however,  costs  of  production 
and  inventory  become  economically  significaunt ,  while 
beyond  this  span,  any  quantitative  approach  rapidly 
approaches  speculation.   The  only  valid  conclusion  then 
would  seem  to  be  that  a  realistic,  quantitative 
approach  to  any  capacity  problem  in  the  production  sec- 
tor must  include  the  effects  of  the  production  and 
inventory  policies. 

A  further  economic  and  realistic  condider- 
ation  in  the  capacity  expansion  problem  must  be  the 
fact  that  capacity,  in  whatever  form  it  may  take,  is 
usually  obtainable  only  in  fixed  increments.   Further- 
more, an  Incremental  acquisition  of  capacity  could,  by 
the  nature  of  the  machine  or  whatever  that  is  added, 
correspond  to  a  large  incremental  capability  in  terms 
of  units  of  output.   Thus,  although  the  roundoff  of  an 
exact,  non-integer  solution  may  give  an  answer  that 
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Is  "cl©se"  to  optimal  in  terms  of  the  demand  profile, 
the  error  in  terms  of  cost  could  be  very  large.   It 
must  be  concluded  then  that,  at  least  in  the  context  of 
capacity  expansion,  integer  solutions  are  important. 

Perhaps  the  classic  early  delineation  of  the 
capacity  expansion  problem  was  presented  by  Marine 
(Manne   6l)  in  his  analysis  of  the  basic  problem  for 
both  deterministic  and  stochastic  demands.   Early 
linear  programming  models  of  the  problem  were  pre- 
sented by  Fetter  (Fetter   6l)  and  Veinott  and  Wagner 
(Veinott-Wagner  62).   An  equivalent  form  of  the 
lattera'  model  and  a  corresponding  network  formulation 
can  easily  be  obtained  from  P2  by  setting  production 
and  inventory  costs  to  zero.   Kalro  and  Arora  (Kalro- 
Arora  70)  extended  the  basic  linear  programming  form- 
ulation to  handle  stochastic  demands. 

A  dynamic  programming  approach  to  capacity 
expansion  is  given  by  Howard  and  Nemhauser  (Howard- 
Nerahauser  68)  for  deterministic  demands  ajid  extended 
to  the  stochastic  demand  case  by  Montgomery  (Montgomery 
71).   An  extension  of  the  dynamic  programming  model 
for  deterministic  demands  to  handle  concave  costs  of 
expansion  is  made  by  Manne  and  Veinott  (Manne-  Veinott 
67),   None  of  the  above  work  takes  into  account  the 
production-inventory  problem,  however,  and  thus  will 
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not  be  considered  as  directly  antecedent  to  this  thesis. 

Zabel  (Zabel  63)  was  one  of  the  early  authors 
to  look  at  the  total  problem,  but  his  assumption  of  no 
inventory  effectively  reduces  the  problem  he  solves  to 
one  of  capacity  expansion.   Hinomoto  (Hlnomoto  65) 
essentially  retained  Zabel's  no-Inventory  assumption 
and  extended  his  work  to  Include  technological  Impro- 
vement, thouch  he  couches  the  work  In  terms  of  replace- 
ment theory  (falling  to  mention  Zabel). 

Iglehart  (Iglehart  65)  seems  to  have  analyzed 
the  first  complete  model  for  the  CPI  problem.   He  ex- 
tends Zabel's  work  to  allow  Inventory  and  restricted 
back-ordering  using  analytical  techniques  to  establish 
optimal  capital  accumulation  and  production  policies. 
He  poses  what  seems  to  be  a  serious  restriction,  how- 
ever, when  he  assumes  that  the  consumptive  factors  of 
production  are  always  spread  evenly  over  all  fixed  cap- 
ital, in  effect  eliminating  the  possibility  of  idle 
capacity.   Carroll  (Carroll  68)  extended  Iglehart's 
work  to  stochastic  demands  and  pricing  considerations 
using  dynamic  and  non-linear  programming  techniques  to 
establish  optimal  capital  Investment  and  production 
policies. 

The  most  significant  general  formulation  of 
the  CPI  problem  to  date,  from  a  quantitative  viewpoint. 
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has  been  the  work  of  Coblan  and  Lele  (Coblan-Lele  70). 
They  give  a  general  dynamic  programming  formulation  of 
a  multi-plant,  multi-Inventory  facility  firm  faced 
with  deterministic  demands  and  show  how  the  model  can 
be  hierarchically  decomposed  to  simplify  the  solution 
process.   Unfortunately,  even  when  only  one  plant  and 
one  Inventory  facility  are  assumed,  the  model  still 
contains  two  state  variables  and  gives  non-Integer  solu- 
tions , 

The  area  of  recent  work  which  appears  most 
relevant  to  the  area  of  capacity  expansion  seems  to 
have  been  done  In  the  related  literature  on  employment 
smoothing.   Interpretatlonally ,  work  force  level  can 
be  viewed  as  a  measure  of  capacity;  the  assumption 
made  In  this  literature  being  that  the  fixed  capacity 
level  Is  given  In  each  period,  as  opposed  to  assuming 
that  the  work  force  Itself  Is  optimally  adjusted  to 
the  level  of  output,  as  mentioned  In  the  previous  sec- 
tion on  economics.   The  greatest  objection  overcome  In 
the  employment  smoothing  literature  Is  that,  by  nature 
of  the  problem,  the  production-Inventory  policy  Is 
always  Included,   To  date,  however.  Integrality  has  not 
been  considered  an  Important  aspect  of  the  problem. 
An  excellent  delineation  of  the  problem  and  methods  of 
solution  plus  a  bibliography  through  1967  Is  given  by 
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Silver  (Silver  67), 

A  linear  programming  model  which  is  close  in 
form  to  P2   and  a  heuristic  for  "rounding-off "  the 
linear  programming  solution  are  given  by  O'Malley, 
Elmaghraby,  and  Jeske  (O'Malley  et  al.  66).   As  is  well 
known,  and  demonstrated  by  example  later  in  this  thesis, 
rounding-off  can  be  far  from  optimal  and  a  very  costly 
approximation  to  the  optimal  integer  solution. 

The  major  difference  between  P2   and  the 
model  of  O'Malley  et  al,  is  the  non-stationarity  of  all 
costs  in  the  former,  as  opposed  to  the  complete  station- 
arity  of  costs  in  the  latter.   In  addition,  the  capa- 
city decrease  coefficient  in  P2   can  be  positive  or 
negative,  within  the  bounds  of  the  preconditions,  to 
correspond  to  salvage  value  or  scrapping  cost  re- 
spectively, while  the  corresponding  term  in  the  O'Malley 
model  is  a  firing  cost  and  constrained  to  be  positive. 
Lastly,   P2  is  an  integer  program  as  opposed  to  the 
linear  program  of  O'Malley.   Still,  O'Malley 's  model  is 
not  a  special  case  of  P2   since  it  has  a  term  for  over- 
time production.   This,  however,  is  a  simple  extension 
of  P2   done  in  Chapter  III,  at  which  time  the  O'Malley 
model  will  become  a  special  case  (except  for  integral- 
ity). 

The  balance  of  the  recent  major  work,  in 
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employment  smoothing  has  been  done  In  two  companion 
papers  by  Llppman,  Rolfe,  Wagner,  and  Yuan  (Llppman  et 
al.  67a  and  67b)  and  a  later  paper  by  Llppman  and 
Yuan  (Llppman-Yuan  69).   The  model  addressed  In  the  1967 
papers  Is  a  generalization  of  the  linear  model  of 
0*Malley  et  al.  with  convex-like,  non-stationary  pro- 
duction costs.  Increasing  and   non-stationary  inventory 
costs,  and  v-shaped,  stationary  smoothing  costs.   De- 
mand-dependent upper  and  lower  bounds  on  work  force 
siie  are  given  find  optimal  policies  are  developed  for 
the  special  cases  of  monotone  increasing  and  monotone 
decreasing  demands.   Algorithms  are  then  given  for  the 
special  cases  of  monotoniclty  and  the  added  restriction 
that  all  costs  are  linear,  though  non-stationarity  is 
retained  for  i>roduction  and  Inventory  costs. 

In  the  1969  paper,  similar  analytical  work 
is  done  on  a  model  which  is  the  same  as  that  of  1967 
with  the  exception  that  smoothing  cost  is  permitted  to 
be  non-stationary.   This  paper  also  points  out  that  the 
upper  bounds  and  algorithms  developed  in  the  1967 
papers  are  invalidated  under  the  more  generalized 
smoothing  costs.   No  further  algorithms  are  presented. 

The  models  in  this  thesis  are  special  cases 
of  the  work  of  Llppman  et  al,  only  Insofar  as  the  more 
general  form  of  the  latter 's  production  and  inventory 


I  •      r         .  f   -  » 


26. 
cost  functions  are  concerned.   Integrality  aind  the 
existence  of  algorithms  In  this  work  seem  to  make 
further  comparison  unnecessary.   The  most  cogent  point 
In  the  above  discussion  would  seem  to  be  the  recogni- 
tion of  the  Interrelationship  between  capacity  expan- 
sion literature  and  employment  smoothing  literature, 
about  which  more  will  be  said  In  the  final  chapter. 

1.6   Plan  of  the  Work 

In  Chapter  II,   PI  and  P2  will  be  trans- 
formed to  the  equivalent  problems  on  which  the  algor- 
ithms given  In  the  chapter  are  based.   Bounds  sire  de- 
veloped here  for  the  variable  which  acts  as  a  surrogate 
for  capacity.   Finally,  more  efficient  bounds  and 
algorithms  are  presented  for  special  monotone  demand 
and  cost  cases. 

Chapter  III  presents  a  network  formulation 
of  the  entire  problem  for  the  special  case  when  cap- 
acity Increase  and  decrease  costs  are  equal  In  each 
period,  but  are  allowed  to  be  non-stationary.   The 
basic  models  are  then  extended  to  Include  Indivisibil- 
ity in  acquisition  of  capacity  and  plecewise-linear, 
convex  production  costs  that  provide  for  consideration 
of  multi-shift  work  and  overtime. 
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A  discussion  of  similar  integer  programming 
algorithms  and  a  reinterpretation  of  the  results  of 
this  thesis  for  the  employment  smoothing  problem  are 
presented  in  Chapter  IV.   Also  presented  here  are  the 
results  of  test  problems  run  with  a  FORTRAN  IV  coding 
of  the  algorithm  along  with  a  summary  and  conclusions. 


Chapter  II 
II. 1  Introduction 

Models  PI  and  P2  of  Chapter  I  could  be  re- 
formulated to  eliminate  the  Inventory  variable  v.  In 
which  case  PI  would  have  2N  constraints  and  2N  variables 
and  P2  would  have  3N  constraints  and  3N  variables.   In 
the  capacity  expansion  context  with  a  period  of  one 
year  and  N  =  5,  neither  problem  would  be  very  large  and 
existing  all  Integer  algorithms  could  be  used.   On  the 
other  hand,  In  the  employment  smoothing  context  with  a 
period  of  one  month  and  N  =  24,  where  furthermore  the 
CIDPI  model  Is  most  relevant,  the  problem  size  Is  out 
of  the  range  of  all  except  branch  and  bound  type  al- 
gorithms. 

Use  of  Implicit  enumeration  seems  out  of  the 
question  since   the  variables  of  PI  and  P2  can  take  on 
rather  large  values,  and  expansion  to  0  -  1  variables 
would  quickly  exceed  the  capacity  of  existing  algorithms 
for  even  small  values  of  N.   Primal  and  dual  cutting- 
plane  algorithms  have  the  disadvantages  of  highly  un- 
predictable convergence  and,  for  the  latter,  lack  of  a 
current  feasible  solution  until  optlmality  is  reached. 
Although  there  is  the  possibility  that  one  of  these 
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other  methods  possesses  an  extremely  efficient  adapta- 
tion for  this  special  problem,  this  thesis  investigated 
the  most  obvious  approach-  branch  and  bound.   In  this 
chapter,  special  branch  and  bound  algorithms  will  be 
developed  for  PI  and  P2  which  exploit  the  structure  in- 
herent in  these  problems. 

In  section  11,2,  a  linear  transformation  is 
done  on  PI  to  obtain  an  equivalent  form  P3  for  which  an 
algorithm  can  be  more  easily  developed.   The  nature  of 
the  algorithm  is  discussed  subsequent  to  the  transform- 
ation.  When  treated  as  a  function  of  certain  of  its 
variables,  the  optimal  value  of  the  objective  function 
of  P3  has  convexity  properties  over  certain  ranges  of 
those  variables  which  are  established  in  section  II. 3 
by  a  series  of  lemmas,  theorems  and  corollaries.   The 
properties  of  the  integer  program  P3  plus  the  results 
from  section  II. 3  lead  to  the  set  of  bounds  and  op- 
timality  tests  given  in  section  11.^  which  are  sub- 
sequently used  in  the  branch  and  bound  algorithm  of 
section  11,5. 

A  similar  development  is  followed  for  the 
CIDPI  model  P2  in  sections  11,7  and  II. 8,   The  linear 
transformation  used  is  the  same  and  the  resulting 
model  P5  has  properties  similar  to  those  of  P3  both  In 
problem  structure  and  objective  function  convexity. 
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The  modifications  to  the  basic  algorithm  of  section 
11,5  needed  for  P5  are  given  at  the  end  of  section  II. 8, 
For  both  transformed  models  there  are  special 
cases  of  demand  and  cost  structures  which  lead  to  major 
simplifications  in  the  algorithms.   These  structures, 
their  properties,  and  the  attendent  modifications  of 
the  algorithms  are  given  in  section  11,6  for  the  CPI 
model  and  section  II. 9  for  the  CIDPI  model. 

II. 2  Transformation  of  PI 

A  set  of  necessary  and  sufficient  conditions 
for  the  feasibility  of  a  w  can  be  obtained  from  con- 
straint sets  (Pl.l)  and  (PI. 2),   The  method  in  general 
is  to  solve  the  inequalities  (Pl.l)  for  x  by  adding 
the  appropriate  terms  to  each  side.   Then  for  any 
k  *  1,2,...,N,  sum  the  first  k  constraints  of  (PI  .2)  to 
obtain  a  constraint  of  the  form 


x^  +  X2  +  ...  +  x^  -  Vj^  =  J^di  -  Vq 


and  substitute  for  the  x.  from  (Pl.l)  giving  an  in- 
equality of  the  form 


>   k 
kw-j^  +  (k-l)w2  ■*■•••■'■  ^ic  ~  "^k  ~   ^  ^i  "•  '^o  ~  ^"i 


-     4         ^    ?JS 


I    •   -     -  ;,  1 
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Repeating  for  each  value  of  k  gives  the  set 


w^  -  v^  -  d^  -  Vq  -  % 


(2.1) 


2wi  +  W2  -  V2  -  d^  +  d2  -  Vq  -  2w^ 


Nw^  +  ...  +  Wfj  ~  ^N  "•  ^1  "^  •••  "^  ^  "  '^o  - 


Nw 


If  a  vector  of  surplus  variables  t  -  0  is 
added  to  (2,1)  and  the  v^  and  t.  in  each  constraint  of 
(2.1)  are  aggregated  as 


Yi  =  v^  +  ^ 


1  =  1,... ,N 


then  y  -  0  since  v  -  0  and  t  -  0,  and  (2.1)  can  be 
written  in  vector  form  as 


(2.2) 
where 


Bw  -  y  =  d  -  Vq  -  w 


B  = 


1 

0 

0 

...00 

2 

1 

0 

...00 

3 

2 

1 

...00 

.   .   • 


N  (N-1)  (N-2)  ...21 


and  the  form  of  the  vectors  is  obvious. 


Solving  for  w  from  (2.2)  gives 
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where 


,-1:: 


w  -  B-^d  +  B~"'-y  -  B-lv„  -  B-lw, 


B' 


-1 


1 

0 

0    0.. 

.    .    0 

0    0 

2 

1 

0    0., 

.     .    0 

0    0 

1 

-2 

10., 

.     .    0 

0    0 

0 

• 

1 

• 

-21., 

.        a        •        < 

.     .    0 

0    0 

•         • 

0  0  ...  1  -2  1 


Expansion  gives  the  following  set  of  relations 


w^  =  d^  +  y^  -  Vq  -  ^o 


(2.3) 


w^  -  d^  -  d,  +  yp  -  2y-  +  v_ 


W3  =  d^  -  d2  +  y3  -  2y2  +  y^ 


\  =  ^  -  ^N-1  +  ^N  -  ^^N-l  "■   yN-2 


which  expresses  w  explicitly  In  terms  of  demand. 
Initial  Inventory,  Initial  capacity,  and  the  newly  In- 
troduced aggregate  variables  y.   The  structure  of  the 
right-hand  side  Is  rather  unique  and,  as  will  be  seen 
shortly,  extremely  useful.   The  above  transformation 
will  be  referred  to  as  Staindard  Transformation  1. 


i'  >  I 
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An  Interesting  economic  interpretation  can  be 
given  to  the  variable  y.   If  overcapacity  is  defined  as 
the  cumulative  excess  of  capacity  over  cumulative  de- 
mand up  to  and  including  period  i,  then  y.  is  the  over- 
capacity in  period  i.   For  instance,  ifN  =  3,   d'^»(3> 
5,6),  and  w*^  *  (5,0,0),  and  w  =  v^  =  0  then  y^  =  y2  =  2 
and  y-^  =  1.   Further,  in  any  period  i,  v.  is  the  por- 
tion of  overcapacity  y^^  which  is  used  to  produce  for 
demand  after  period  i  and  stocked  as  inventory,  and  t^ 
is  the  cumulative  amount  of  idle  overcapacity  as  of  the 
end  of  period  1,   Thus,  if  Vj.  »  0,  then  yj^  is  the  total 
idle  capacity  in  the  model. 

Model  PI  can  be  transformed  by  substituting 
for  w^  from  Standard  Transformation  1  to  give  the 
following  equivalent  model, 

P3:     min  Z^  =      l   Lp^x^   +  h^v^^  +  c*yj^]  +  K 


3.  t. 


(P3.1) 


X,  -  dT  +  y.  -  V, 


^i  "  ^i  "*■  ^i  '  ^i-l     ^  *  2,...,N 


■'•For  the  remainder  of  this  thesis,  the  objective 
function  of  ?^   will  be  referred  to  as  Z^, 
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^  -  ^  ■  ^1  -  ^o 


(P3.2)         x^  +  Vj^_^  -  Vi  «  dj,      i  =  2,...,N-1 


*N  *  ^N-1  =  ^N  +  Vj^ 


(P3.3)    72  -  2y^  "  ^1  ~  ^2  -  '^o 

^i  ~  2^1-1  "^  ^1-2  ~  ^i-l  "  ^i     i  =  3,...,N 
3c,  V,  y  -  0    and  INTEGER 


where  constraint  set  (P3.3)  comes  from  substitution  in 
w  -  0  in  PI;  the  restriction  of  y  to  integers  comes  from 
the  restriction  of  w  to  integers  in  PI  and  the  fact 
that  the  right-hand  side  of  constraint  set  CP3.3)  is  in- 
teger by  assumption,  and  where 


c^  =  Cj^  -  2cj^^^  +  Cj^+2  +  ra^  -  n^i+i   1  "  l,..,,N-2 
jj_l  "  ^N-1  -  2cjj  +  n^_-L  -  m^ 


c 


N-1 
K  -   ^^    [(c^-Ci^.^+mj^).dj^]  +  (cj^+mj^)djj-C3_WQ+(c2-ci-m3^)vQ 


If,   in  P3,    a   feasible  y   is    given   then  the  problem 


^      f 
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becomes 


N 
?H:  mln   Z^   =      Z    Ip^^:^   +   ^l^i^ 


s.  t. 


CP^-1) 


^1  -   ^1  ^  ^1  -  ^o 


^1  ~  ^1  ■•■  ^1  -  yi-i       ^  "  2,...,N 


1  1  1  o 


(P4.2)  Xj^   +   Vj^^^  ~   ^1   '   ^1  ^   '   2,...,N-1 


X,    V  -  0        and  INTEGER 


which  Just   is   the  production-inventory   problem,   para- 
metric  in  y;    but   which,   more   importantly,    can  be  shown 
to  be   a  tramsportation  problem  as   was    first  pointed  out 
by  Bowman    (Bowman   56).      Specifically,    let 

^il"  production   in  period  i   for  demand 
in  period  J,   J    -  i 
then   P^   can  be  written   as 


N        N 
?H*  :         min  z;    =     Z        2   a.  .xii 
4        i=l  j«i   iJ    ^^ 


',■■    •<• 


.        >      .     ^ 
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3.  t. 

CP^.l')      ?  XII  -  J   "•        ^ 

J-i       [d^   +  Yi  -  71-1      i  "  2,...,N 


J 


/" 


di  -  Vq  J  ■  1 


(P^.2')        l  xi^  -  '  dj  J  -  2,...,N-1 

i"!       dM  +  VxT  .1  -  N 


l^clN  +  Vjj  J  -  N 


X  -  0    and  INTEGER 


where 


J-1 
aij  »  Pi  +  ^l^h^  J  -  1+1,.... N 

ail  ■  Pi 


and  X  Is  Integer  by  the  properties  of  the  transporta- 
tion problem. 2 

The  impact  of  this  development  Is  that.  In  a 
branch  and  bound  algorithm  enumeration  would  have  to 
be  done  on  only  a  subset  of  the  problem  variables,  i.e. 
y.   For  each  feasible  y  generated,  corresponding  to  a 
node  in  the  enumeration,  a  transportation  problem  can 
be  solved.   This,  in  effect.  Implicitly  fathoms  each  time 
what  would  have  been  a  significant  portion  of  the 

^The  lower  summation  range  in  (P4.1M  could  be  ex- 
tended to  1  and  the  upper  range  in  (PA.2')  to  N  if 
aj^j  =00,VJ  <  i,  1  =  2,...,N. 
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enumeration  tree  had  branch  and  bound  been  tried  on 
the  entire  problem.   Furthermore,  the  number  of  trans- 
portation problems  solved  can  be  cut  to  a  small  fraction 
of  the  number  of  nodes  actually  visited  by  using  suit- 
able tests  and  bounds.   In  effect,  the  branch  and  bound 
is  a  search  for  the  minimum  cost  combination  of  a 
feasible  y  and  its  corresponding  optimal  production- 
inventory  schedule. 

Lower  bounds  for  y  can  be  obtained  directly 
from  constraint  set  (P3.3)  which  can  be  solved  for  y  as 

^1  -  «o  -^  ^o  -  ^1 
(2.4)     y2  -  2y^  +  d^  -  d^  -  v^ 

y^  -  2yj^_^  -  yj^_2  +  '^i-l  -  ^i     i  =  3,...,N 

This  set  of  lower  bounds  on  y  is  integer  since  all 
terms  on  the  right-hand  side  are  integer. 

Upper  bounds  on  the  problem  variables  must 
exist  since  the  model  is  bounded  .   The  following  lemma 
will  give  a  set  of  global  upper  bounds  on  y. 
Lemma  2.1  Global  upper  bounds  for  y  are 


<    N       J 
(2.5)     y.  -  J-  j:  d.  -  Ed.-  (J-l)v.    J  -  1,...,N 
J     i«l  ^   i»l  ^         ° 


1   ■■    '.".': ...^  J  -.  •e  '.  ''^  •    H  '.' 


'•i:         '  r 


'  N      - 
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Proof:   For  J  =  1  from  (2.2) 

y^  -  w^  -  d-L  +  Vq  +  Wq 

and  by  Inspection  the  right-hand  side  Is  maximized  when 
w-j^  takes  on  its  maximum  value.   The  maximum  capacity  In 
period  1  would  be  a  level  sufficient  to  produce  all 
future  demands  In  period  1.   This  value  Is 


N 
and  so 


Wt  =   L  d^  -  v^  -  w^ 
1   jL=l  1    o    o 


max 
or 


^1  =  Ja  -  ^1 


N 


^1  -  i=i^i  -  ^1 

In  general  let  J  =  k,  ke{2 ,3 , , . . ,N} 

k 
^k  ~  ^^1  "*"  ^^""1)^2  +  . .  .  +  wi^  -   ^  ^1  +  Vq  +  kWQ 

Again  the  right-hand  side  Is  maximized  when  w^^  Is,  and 
from  the  argument  above 

<    N       k 

y,  ^  k.  E  di  -   Id.  -  (k-l)v 
^     1»1     1»1  ^         ° 

so  the  result  Is  true  In  general. 

QED 


.  k.  :•.    ,.    V  o  ^    (, 


.J' 


-  -J  , 
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As  will  be  discussed   later  in   this   chapter, 
upper  bounds    are  not   really    crucial   to   the   algorithm 
since  more   powerful   termination  rules  will  be   employed. 
A  special   case  which  has   an   optimal   solution  with 
yj^   -   0  will  also  be   investigated  subsequently. 

II. 3    Properties   of   P3 

Since   a  branch  and  bound  algorithm  is   to  be 
developed  it   is  important   to  have   efficient  bounds   and 
termination  rules   in   order  to  perform  as    little   enumer- 
ation as  possible.      The   following  lemmas  will  establish 
useful  properties   of   some  variables   of  P3. 
Lemma  2.2     If  all  y.   are  fixed  for  1  =  1,...,N-1  and  y^  is 
incronented  by  n  to 

then  the  new  value  of  w„  Is 

t/lt'   =  w     +  n 
^N         N 

Proof:     Prom  (2.3) 


=   ^V"^   -  2%-l  ""  ^N-2  ^  ^^  -  Vl 


=  ^N  "^  "^  QED 


For  Lemmas  2.3  and  2.4,  Corollary  2.2.2  and  Theorems 
2.1  and  2.2,  when  given  a  set  of  y. ,  1  =  1,...,J,  It  Is 


^0, 


assumed  that  the  remaining  y         are  at  their  lower  bounds  for 
k  =  1,...,N-J. 


Lemma   2,3      If  for   given  y . ,    i    »    1, . . . ,j ,    the    lower 
bounds    of  yj+k   ^^^  positive,    k  =    1,...,N-J    then  wj+j^   *   0, 
k   -    1,. . . ,N-J . 

Proof:      For  any   ke{l,2, . . . ,N-J }    from   (2.^)    and   the 
hypothesis 

yj+k  -   2yj+it.i  -  yj+k-2  +   ^i+k^l  -   «^J+k   ^    ^ 

Further,    if  yj+k  i-s   ^^   i^s    lower  bound,    equality   holds 
on   the   left.      The   expression   for  w.^,     from   (2.3)    is 

''j+k  "   yj+k  -    2yj+k-l  +   yj+k-2  +    dj+k-1   -   dj+k 

but   the   right-hand  side  has   Just   been   shown   to   be   zero. 
Thus 

QED 
Lemma   2 « U      If  for   given  y    ,   i   «    1,...,J,    the    lower 
bounds   of  yj+k   si^®  positive,    k  »    1,...,N-J,   and  yj    » 
y.   +   n,   ne{l,2,3, . , , }    then 


w*    «  W4   +  n 

1  J 
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and  for  k  «    1,2,,..,N-J,   the   lower  bound  on  y'        is 

yj+k  +    (k+l)n 


and 

w'        =»   0 

Proof:      From   (2.3) 

Wj    =  yj    -   2yj_i  +   yj_2  +   d^    -   dj_i 


so 


wj    -    (yj+n)    -   2yj_^  +  yj_2  +   ^j   -   dj.^  »  Wj    +  n 

which  proves  the  first  assertion.   The  second  and  third 
are  proved  by  Induction  on  k.   For  k  =  1,  from  (2.4) 


and  from  (2,3)  and  Lemma  2.3 

Assume  the  relationship  holds  for  a  general  ke{2,3,..., 
N-J-1}.   Then  for  k+1  from  (2.4) 
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=  yj+k+1  +  (k+2)n 

The  final  assertion  Is  then  true  by  Lemma  2,3. 

QED 
Remark:   From  the  relationships  (2.3)  and  (2.U)  It  can 
be  seen  that  a  change  In  y.  leaves  wj^, , , .  ,Wj__j|^  and  y-,, 
•••fyj-l  unaffected. 
Now  let 

CCy^  «  capacity  cost  as  a  function  of  y. 
PKy^)  "  production-Inventory  cost  as  a 
function  of  y^ 
T(y.)  *  total  cost  as  a  function  of  y. 
so 

T(yj)  «  C(yj)  +  Pl(yj) 

Then  the  preceedlng  lemmas  can  be  used  to  establish 
the  following  cost  properties  for  P3: 
Theorem  2. 1   If  for  given  y.,  1  =  1,.,.,J,  the  lower 
bounds  for  y^^^.  are  positive  k  =  1,...,N-J  then  CCy^+n) 
is  linear  and  non-decreasing  in  n  for  ne{l,2,3 . . • }. 

Proof:   From  Lemma  2.4  for  y'  =  y.+n,  n6{l,2,3...} 


w'  «  w.  +  n 
wLjj.  =  0  k  «  1,2,. . .  ,N-J 


1  ^  V ;  i 


I      ;■      o 


n'i 


^3. 
Thus,  from  PI 

C(y  +1)  -  C(y.)  »  C(y  +2)  -  C(y,  +  1)  =...  =  (c.+  E  m.  )  -  0 
J         J       J         J  J  1=J 

QED 
Before  stating  the  next  theorem,  a  well 
knovm  result  for  the  linear  progrsim  associated  with 
P^*  must  be  mentioned.   Since  the  linear  program 
satisfies  the  sufficient  conditions  for  unlmodularity , 
then  all  extreme  points  of  the  convex  hull  are  Integer 
(Hu  69).   This  allows  use  of  the  properties  of  the 
general  linear  program  in  the  following  proof. 
Theorem  2.2   If  ne{l,2,3, . . . }  and  for  any  Je{1,2, . . .  ,N} 
if  yj  =  y*  +  n  and  for  each  value  of  n  the  y.  are  set 
at  their  lower  bound,  1  =  J+1,.,.,N,  then  Pl(y')  is 
convex  in  n  and  the  y . ,  1  =  J+1,.,.,N. 

Proof:   Assume  the  following  linear  programming 
problem 

(pi):        min  c^x 
s,  t. 

Ax  =  b 


X  i  0 


and  the  two  pararaeterizations  of  b 


^^1 
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(p2)  Ax  -  b  +  (l-X)e 

(p3)  Ax  =  b  +  e 

where  0  -  X-  l,e  is  any  vector  which  maintains  feasi- 
bility in  (p3)  and  the  optimal  solution  of  (pi)  is  x^^ 
and  of  (p3)  is  x^  where 

c'^x-L  =  f(b) 

c'^Xo  »  f(b+e) 

are  the  optimal  objective  function  values  as  functions 
of  the  right-hand  sides.   Now  let 

(2.6)  x»  =  Xx,  +  (l-X)x^     0  -  X  -  1 

and  since  the  problems  are  linear  it  must  be  that 


(2.7)  c'^x^   =  c"^^  +  (l-X)c'^x 


Multiplying  the  constraints  of  (pi)  and  (p3)  by  the 
appropriate  constants  gives 


XAi.  =  Xb  ^   ^ 

^  0  -  X  -  1 

(l-X)Ax   -  (1-X)(b+^) 


^5. 
and  since  X  Is  a  scalar,  this  can  be  rewritten  as 


AXxt  =  Xb 

0  -  X  -  1 


AC1-X)X2  »  (1-X)(b+^) 


Adding  these  two  constraint  sets  gives 

A[XXj_  +  (1-X)X2]  =  b  +  (l-X)e 
which,  from  (2.6),  Is  Just 


Ax'  -  b  +  (l-X)e 


which  proves  that  the  convex  combination  of  the  op- 
timal solutions  to  (pi)  and  (p3)  Is  feasible  to  the 
convex  combination  of  right-hand  sides.   Now  if  x  Is 
the  optimal  solution  to  (p2) 


(2.8)  c^x^  -  c'^x^ 


since  It  Is  a  minimization  problem. 


From  (2.7)  and  (2.8) 


c^X2  -  Xc'^i  +  (l-X)c'^x-j     0  1  X  -  1 


-.'  J  ^ 


'   )      -    0     » 


'/_  --•  <•    I  - 
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producing  the  general  convexity  result  and,  by  the 

remark  preceding  this  theorem,  the  desired  result  for 

PI(yj+n). 

QED 

The  following  corollaries  follow  trivially. 

Corollary  2.2.1   PKy^^+n)  is  convex  for  n£{l,2,3, . . .  )  . 

Corollary  2.2.2   If  for  given  y^ ,  1  =  1 , . . . , J ,  the 

lower  bounds  of  yj+j^  are  positive,  k  =  1,...,N-J  then 

PlCy,*n)   is  convex  in  n  for  ne{l,2,3, . . . } . 

Proof:   From  Lemma  2,H   the  capacities  In  periods' 

J+1  through  N  are  equal  to  the  capacity  in  period  J. 

The  result  follows  from  this  and  Theorem  2.2.   QED 

The  major  result  can  now  be  given  as  follows; 

Theorem  2.3   If  for  given  y.,  1  =  1,...,J,  the  lower 

bounds  for  y^+j,  are  positive,  k  =  1,...,N-J,  then 

T(y.+n)  is  convex  in  n  for  ne{l,2,3, . . . } . 
J 

Proof:   The  proof  follows  from  Theorem  2.1  and 
Corollary  2.2.2  and  the  fact  that  the  sum  of  convex 
functions  is  convex.  qed 

Corollary  2.3.X  T(yj^+n)  is  convex  for  n£{l,2  ,3, .  . .  } . 

II. 4   Qptimality  Tests  and  Bounds  for  the  Enumeration 

The  results  of  section  II. 3  can  now  be  used 
to  develop  bounds  and  optimallty  tests  for  the 


algorithm.   The  logical  first  step  is  to  solve  the 
linear  programming  problem  for  PI  and  use 
Optimality  Test  1_:   If  the  optimal  linear  programming 
solution  is  all  Integer,  then  this  is  an  optimal  solu- 
tion to  P3. 

Otherwise  let 

Z*   =  optimal  objective  function  value  of  the 
LP 

linear  program  associated  with  the  integer 
programming  problem 
Z*   =  optimal  integer  program  objective  function 
value 
It  is  a  well  known  fact  that 


Z»  ^   Z» 

LP     IP 


and  this  is  used  as 

Bound  1:   A  lower  bound  on  total  cost  is  Z*  . 
LP 

The  algorithm  will  use  Bound  1  in  the  follow- 
ing optimality  test. 

Optimality  Test  2:   If  Zjp  is  the  current  integer  pro- 
gramming solution  and 

then  Zjp  is  an  optimal  solution  to  the  integer 


.  A   I'l.      ^    .-  C 


•M 
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program,^ 

The  next  bound  uses  the  minimum  production- 
inventory  cost,  PI*,  to  obtain  an  upper  bound  on 
capacity  expansion  cost.   PI»  can  be  obtained  by 
solving  P4'  with  the  values  of  the  right-hand  sides 
of  constraint  set  (P^,l')  equal  to  infinity.   This 
corresponds  to  putting  an  infinite  supply  at  the 
sources  of  the  transportation  problem. 

Lemma  2,8   If  Z*  is  a  current  optimal  solution  in  the 
3 

enumeration  then 

a  =  Z'  -  PI» 
3 

is  an  upper  bound  on  capacity  expansion  cost. 

Proof:   Assume  there  is  a  capacity  expansion  cost 


a   such  that 


1   " 
a-^  >  a 


and  with  a  corresponding  Z^   such  that 

3 


.1      .      rr, 


z-^  <  z 
3         3 


But   this   implies   that 


3 Pa]    indicates  the  next    integer   greater   than  a, 


^9 


PI  «  Z"^  -  a^  <  Z'  -  a  =  PI» 


which  is  a  contradiction  of  the  construction  of  PI*. 

QED 

Using  the  above  gives: 

Bound  2:   If  Z *  is  a  current  optimal  solution  of  the 
3 

integer  program,  then  an  upper  bound  on  capacity  ex- 
pansion cost  a  is 


a  -  Z*  -  PI«  =  a 
3 


Lemma  2,8  can  be  useful  in  another  way.   If 
the  current  upper  bound  on  capacity  cost  is  a  then 
certainly  an  upper  bound  on  period  1  capacity  is 


< 

"l  - 


N 
Ct+  L   m^ 

-L  1=1  1 


since  any  larger  value  of  w-j^  would  violate  a.   A 
similar  bound  for  W2  is 


W2  - 


a-Cci +  Z  m. )w, 
J-  1=1  1   1 

N 


laj  indicates  the  greatest  Integer  smaller  than  a 
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In  general  the  following  bound  can  be  used. 
Bound  3:  For  the  current  a,  upper  bounds  on  capacity 
expansion  In  each  period  are 


< 

^1  - 


< 

"J  - 


a 


N 
Ci  +  Z   m^ 


1=1   ^  k=l  ^   ^ 

N 
c  4  +  I  mi. 
J  k«J  ^ 


J  »  2,...,N 


The  following  bounds  use  the  results  of  the 
previous  sections. 

Bound  k_:      For  any  feasible  values  of  y.  ,  1  »  1,...,J, 
lower  bounds  on  the  values  of  y.,  1  «  J+1,...,N  can  be 
obtained  from  (2.4). 

Bound  5_:   Global  upper  bounds  on  y.,  1  «  1,...,N  can 
be  obtained  from  (2.5). 

Qptlmallty  Test  3^:   If  enumeration  has  been  completed 
with  y-,  at  Its.  global  upper  bound  or  the  upper  bound 
Inferred  from  Bound  3,  whichever  Is  lower,  then  the 
current  optimal  solution  Is  optimal  for  the  Integer 
program. 

The  following  bound  on  yj^  uses  the  results 
of  Corollary  2.3.1. 
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Bound  6_:      If  for  a  given  set  of  y .  ,  1  =  1,...,N-1 

then  the  enumeration  on  y»,  for  this  set  of  y^^  can  be 
fathomed  at  (y^r+i;. 

The  following  bound  uses  the  result  of 
Theorem  ^.3. 

Bound  ]_:      If  for  a  given  y^,  y2*'«*»yj»  ^^^  lower 
bounds  on  y^,  are  positive  for  1  =  J+1,...,N  and  If 

T(yj)  -  T(yj+1) 

then  the  enumeration  on  y^  can  be  fathomed  for  this  set 
of  y^^,  1  =  l,...,j-l  at  (yj  +  i;. 

Bounds  8  and  9  make  use  of  Theorem  2.1  and 
Corollary  2.1.1. 

Bound  8:   For  given  y^^,  1  =  1,...,J,  If  the  lower  bounds 
on  y^^  are  positive  for  1  =  J+1,...,N  and  If  C(y^)  -  a, 
then  y^    Is  fathomed  for  the  preceding  set  of  y^^. 
Bound  9:   If  for  given  y^^,  1  =  1,...,N-1,  CCyjj)  -  a, 
then  no  further  enumeration  on  y»T  Is  necessary. 

Since  the  algorithm  conducts  the  search  by 
enumeration  on  yjj  most  often.  Bounds  6  and  9  become  very 
powerful  tools.   Bounds  7  and  a   are  also  extremely 
powerful  because  their  conditions  will  be  satisfied 


i:  Li'       ^f 
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often  during  enumeration  on  the  later  period  variables 
due  to  relationships  (2,4),   The  usefulness  of  these 
bounds  diminishes  the  need  for  the  global  upper  bounds 
(2.5)  since  such  bounds  must  Include  all  cases  where 
Bounds  6  through  9  would  be  useful, 

II. 5  Algorithm  for  P3 

Step  £  Initialization  Set  the  value  of  current  optimal 
cost,  Z',  and  the  current  bound  on  capacity  cost,  a,  to 
Infinity. 

Step  1   Solve  the  associated  linear  program.   If  the 
solution  Is  Integer,  stop,   (Optlmallty  Test  1)   Other- 
wise set  the  linear  programming  bound  to  (Bound  1) 


P-'1 


and  go  to  Step  2a. 

Step   2a      Solve  an  uncapacltated   transportation  problem 

to    get    the  minimum  production-Inventory   cost    PI*   and    go 

to    Step   2b, 

Step  2b  Find  global  upper  bounds  for  y^,  i  =  l,,,,,N 

from  (2,5).   Oo  to  Step  3. 

Step  3_  Start  the  enumeration  by  setting 
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y^  =  max  [ (w^+v^-d^ ) , 0] 

and  finding  lower  bounds  for  the  remaining  values  of  y^^, 
1  =«  2,...,N,  from  (2.4),  using  this  value  of  y^ .   The 
pointer  (PNT)  is  N.   Go  to  Step  5. 
Step  ^  PNT  =  PNT  -  1.   Go  to  Step  Mb. 

Step  4b  With  PNT  =  J,  let  y^  =  y.+l  and,  If  J  <  N,  find 
new  lower  bounds  for  succeeding  y^  from  (2.4).  If  J  /  1 
and  any  lower  bound  exceeds  an  upper  bound  (from  Bounds 
3  and  5),  yj  is  fathomed  so  go  to  Step  4a  with  PNT  =  J. 
If  ,j  =  1  and  some  lower  bound  exceeds  an  upper  bound,  go 
to  Step  10  (Optimality  Test  3).  Otherwise  go  to  Step  5 
with  PNT  =  J. 

Step  5   With  PNT  =  J  determine  capacity  cost  a  =  C(y.). 
If  a  >  a  and  either  j  «  N,  or  J  <  N  and  y^^  >  0,  1  *  J +1 , 
...,N  go  to  Step  4a  with  PNT  =  J  (Bounds  8  and  9). 
Otherwise  if  a  >  a  go  to  Step  4b  with  PNT  »  N.   If 
a  ^  a  and  if  J  <  N  with  y^  >  0,  1  =  J+1,...,N  go  to 
Step  6  with  PNT  =  J  (This  sets  up  convexity  test  for 
Bound  7).   Otherwise  go  to  Step  6  with  PNT  =  N  (This 
sets  up  convexity  test  for  Bound  6). 
Step  6   With  PNT  =  J  solve  a  transportation  problem 
to  obtain  the  related  optimal  production-inventory  cost 
Pl(yj).   Then  total  cost  is  T(yj)  »  C(yj)  +  Pl(yj).   Go 
to  Step  7  with  PNT  =  J. 


r  i. 
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Step   7      With   PNT   «J,    lfJ=NorJ<N  and   y^^    >    0, 

1  »  J+1,...,N  and  TCy.)  >  TCy.-l)  then  y.  is  fathomed 

so  go  to  Step  4a  with  PNT  =  J  (Bounds  6  and  7).   If 

T(y.)  =  T(y.-l)  »  Z'  in  any  case  go  to  Step  9.   For  all 
J       J       3 

other  cases  If  T(y^)  <  Z '  go  to  Step  8,  otherwise  go  to 

Step  4b  with  PNT  =  N. 

Step  8   If  T(yj)  =  \^*\      SO  to  Step  10  (Optlmallty 

Test  2).   Otherwise  go  to  Step  9. 

Step  9   Either  Include  T(yj )  =  Z'  on  the  list  of 

current  optimal  solutions  or  replace  the  current  list 

with  Z'  =  T(y^),  whichever  Is  appropriate.   Update  In- 

3      '^ 
termedlate  bounds  from  Bound  3  and  go  to  Step  4b  with 

PNT  =  N. 

Step  10  Decode  the  current  optimal  solutions  and  STOP. 

The  algorithm,  as  written,  will  also  find 
all  alternative  optimal  solutions  If  termination  Is  not 
due  to  Step  8.   The  flexibility  of  choice  among  solu- 
tions would  normally  be  desirable  In  a  management 
situation.   If  the  alternate  optima  and  longer  running 
times  are  not  desired  then  the  following  changes  to  the 
algorithm  are  needed: 

Step  5.  With  PNT  =  J  determine  capacity  cost  a    =  C(yj), 
If  a  ^  S  and  either  J  =  N  or  J  <  N  and  y^  >  0,  1  =  J+1, 
...,N,  go  to  Step  4a  with  PNT  =  J.   Otherwise  If  a  -  3 
go  to  Step  4b  with  PNT  =  N.   If  a  <  a  and  If  J  <  N  with 
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y^  >  0,  1  =  J+1,...,N  go  to  Step  6  with  PNT  =  J. 

Otherwise  go  to  Step  6  with  PNT  =  N. 

Step  7   With  PNT  -  J,  If  J  =  N  or  J  <  N  and  yj^  >  0, 

1  «  J+1,...,N  and  T(yj)  -  T(yj-l)  then  yj  Is  fathomed 

so  go  to  Step  Ha   with  PNT  =  J.   If  T(y.)  <  Z'  go  to 

•^      3 
Step  8,  otherwise  go  to  Step  kb   with  PNT  =  N. 

Step  2  Replace  the  current  optimal  solution  with 

Z'    =  T(yj).      Go    to    Step   4b  with   PNT   =   N. 

II. 6     The   CPI  Model   Under   Special   Demand  and  Cost 
Assumptions 

In  this   section   two    special   cases  of   the  CPI 
model   will   be  discussed,    one   Involving  monotone  non- 
decreasing  demands  and   the  other  monotone  non-lncreas- 
Ing  demands.      The  former  leads   to    significant   simpli- 
fications  In  the   CPI  algorithm,    while   the  latter   Is  a 
trivial  case  of   the   CPI  model. 

The  following   two   cost   assumptions   hold   In 
this    section: 

(a)  The  cost   assumptions  of   PI   hold. 

(b)  The  cost  coefficients  of   PI    (Chapter   I)   are  mono- 
tone non-increasing  with  time. 

Assumption   (b)   allows   the  possibility  of  discounting. 


I  .      ■;:'.> 
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Monotone   non-decreasing   demand 

The  case  of  demand  which   is  monotone   non- 
decreasing   is   likely    to   occur  often  in  the  context  of 
firms    interested   in   solving  capacity   expansion  problems 
The  main  result   for   this  case   is  contained   in  the 
following   theorem. 

Theorem   2.4      If  assumptions    (a)   and    (b)    hold,    demand 
is  monotone   non-decreasing,    and 


<      N 
Nw^   +  v_   -      Z  d. 
o  o        i=l    1 


then  there   is  an  optimal   solution  with  Yj^   =   0 
Proof:      Assume   in  contradiction   that 


y„  =  a  >   0 


It  can  also  be  assumed,  as  discussed  in  Chapter  I,  that 

v„  *  0.   But  this  Just  Implies  that  there  are  a  >  0 

idle  units  of  capacity  in  the  model.  Consider  the  first 

unit  of  idle  capacity.   One  of  two  cases  can  be  assumed 

to  exist. 

Case  I-  The  idle  unit  was  added  in  period  J  Instead  of 

k  where  it  was  first  used,  where  J  <  k  ^  N. 

Case  II-  The  capacity  was  added  in  some  period  previous 
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to    .^ ,    was  used   up    to   and    including  J-1,    and   became 
idle   in  J . 

For  Case   I   to    happen  either   the  cost   of  adding  capacity 
in  J   must   be  cheaper   than  k,    i.e. 

N  N 

c,    +      I   m^    <    Ci,   +      li   m^ 
J         l=j    1  ^        i=k   1 

which,  since  c^  -  Cj,  (assumption  (b))  Implies  that 

k-1  "^ 

I  m^^  <  0  contradicting  assumption  (a);  or  the  cost  is 

equal  to  that  in  k  in  which  case  the  idle  unit  could  be 
removed  in  J  and  added  in  k  with  no  change  in  cost. 
For  Case  II  it  must  be  that  w^  =  0  since  if  not,  then 
Case  I  applies.   So  this  unit  of  capacity  must  have 
been  built  in  some  previous  period  q  and  used  to  pro- 
duce from  q  through  J-1.   Now  at  least  one  of  these 
(J-Q)  units  of  production  must  be  inventoried  into 
period  J  since  if  vj_;j^  =  0  then  since  x.^-^    >   x^  this 
implies  that  d^  -j.  ^  ^^    contradicting  the  demand  assump- 
tion.  Furthermore,  by  the  same  argument,  v^  -  1,  i  = 
q,q+l, . . . , j-1  so  the  idle  unit  in  J  was  used  to  produce 
for  inventory  in  period  q  and  at  least  that  unit  of  in- 
ventory was  carried  into  period  J .   The  production 
schedule  of  periods  q+1  to  J-1  will  not  change  if  the 
unit  of  capacity  is  added  in  q+1  instead  of  q,  however 
the  unit  which  was  inventoried  from  q  to  J  could  be 
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produced    In  J    instead,    eliminating   the    idle  unit   of 
capacity.      With   idle  capacity  the   incremental   cost  con- 
sidering  .just    the  one  unit   of   production   in  question 
was 


(2.9)  Cq   +   J  rai   +   pq   +   'l\ 


while  adding  the  capacity  in  q+1  and  producing  in  J 
gives 

(2.10)  Oq+i  +  ^J^^mi  *   Pj 

Subtracting  C2.10)  from  (2.9)  gives 


tCq-Cq+l)  +  niq  +  ^Pq-Pj  )  +  .5  ^i  "  ^ 


by  assumptions  (a)  and  (b).   Thus,  in  both  cases,  idle 
capacity  can  be  reduced  with,  at  worst,  no  reduction  in 
cost. 

If  idle  capacity  still  exists,  the  argument 
can  be  repeated. 

QED 
Corollary  2.^.1;   If  the  assumptions  of  Theorem  2. A 
hold  then  all  capacity  is  used  in  each  period. 

The  use  of  Corollary  2.^.1  leads  to  an  imme- 
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dlate  simplification  of  the  production-inventory  prob- 
lem P^ .   Since  the  right-hand  sides  of  constraint  set 
(P4.1)  are  Just  capacities,  equality  holds  in  all  N 
constraints  and  the  x^  are  determined  trivially  for  a 
given  7,   Constraint  set  (P4.ii)  becomes  N  equations  in 
N  unknowns  (since  v   is  explicit  and  v„  is  either  zero 
implicitly  or  set  explicitly)  and  is  thus  deterministic, 
So  the  subproblem  in  the  enumeration  can  be  solved  by 
inspection  each  time. 

Another  observation  follows  from  Theorem  2.4 
and  Lemma  2.4.1,   If  all  capacity  is  always  used  it  can 
certainly  never  be  feasible  to  build  capacity  in  any 
period  sufficient  to  produce  all  remaining  demand.   In 
fact,  the  upper  bounds  on  capacity  increase  are  much 
smaller  than  those  given  by  (2.5).   Instead  of  stating 
them  explicitly,  they  can  be  handled  implicitly  in  the 
algorithm  by  the  following: 
Bound  10;   For  given  y^^,  i  =  1,...,J,  if  the  increment 


yj  =  yj  ^  1 


causes  the  succeeding  lower  bounds  to  contain  yj^  >  0, 
then  yj  is  fathomed  for  the  preceding  set  of  y^^. 

Bounds  6  through  9  are  now  useless,  since 
their  conditions  are  never  met  in  this  case,  so 
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Bound  10  takes  over  as  a  main  fathoming  rule  in  the 
algorithm.   Furthermore,  the  global  upper  bounds  from 
(2.5)  are  Irrelevant  to  this  special  case  and  can  be 
ignored , 

The  basic  algorithm  requires  the  following 
modifications: 

1)  Drop  Step  2b. 

2)  Step  4b'   With  PNT  =  J  let  yj  =  y.+l  and  find  new 
lower  bounds  for  the  succeeding  elements  from  (2.4). 
If  Yj,  >  0  then  y.  is  fathomed  so  go  to  Step  4a  with 
PNT  =  J  (Bound  10),   If  J  =  1  and  y^  exceeds  its  upper 
bound  from  Bound  3  or  y»,  >  0,  go  to  Step  10.   Otherwise 
go  to  Step  5'  with  PNT  =  J. 

3)  Step  5*   With  PNT  =  J  determine  capacity  cost  a  = 
C(y.).   If  a  >  a  go  to  Step  4b'  with  PNT  =  N-1 .   Other- 
wise go  to  Step  6'  with  PNT  =  J. 

4)  Step  6*   With  PNT  =  J,  determine  the  related  optimal 
production-inventory  schedule  by  inspection  and  the  re- 
lated cost  Pl(yj).   Total  cost  is  T(y.)  =  C(yj)  +  Pl(yj). 
Go  to  Step  7*  with  PNT  =  J. 

5)  Step  72.  With  PNT  =  J  if  T(yj)  >  Z'  go  to  Step  4b' 
with  PNT  =  N-1.   If  T(yj)  =  z»  go  to  Step  9'.   Other- 
wise go  to  Step  8. 

6)  Step  9'   Either  include  T(y4)  =  Z '  on  the  list  of 
J     3 

current  optimal  solutions  or  replace  the  current  list 
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with  Z'  =  T(yj),  whichever  Is  appropriate.   Go  to  Step 
4b'  with  PNT  =  N-1. 


Monotone  non- Increasing  demand 

Assume  that  demand  Is  monotone  non-lncreaslng 
and  that  monotoniclty  Is  not  destroyed  by  Initial  In- 
ventory, I.e. 


di  -  V„  i  dj 


By  assumption  (b)  production  should  always  take  place 
as  late  as  possible.   Then  for  the  CPI  model  the  only 
optimal  capacity  expansion  policy  can  be  to  expand 
capacity  In  the  first  period  to  satisfy  the  demand  In 
that  period.  I.e. 


^1  =  ^1  -  ^o  -  ^o 


This  will  then  be  sufficient  to  produce  the  remaining 
demands,  so  no  capacity  expansion  problem  need  be 
solved.   If  the  cost  assumptions  do  not  hold  however, 
then  this  case  Is  Just  handled  with  the  general  CPI 
model , 
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II. 7   Transformation  for  P2 

A  similar  procedure  to  that  used  In  trans- 
forming the  CPI  model  PI  will  be  used  to  transform  the 
CIDPI  model  P2 .   First  a  lemma  will  be  proved  which 
will  be  useful  In  Interpreting  the  transformed  model, 
and  In  the  algorithm  to  follow. 
Lemma  2 . 9  For  any  period  1  either 

(1)  Wj_  >  0  and  r^^  =  0 


or 


or 


(2)  Wj^  =  0  and  r^^  =  0 


C3)  w^  =  0  and  r^^  >  0 
Proof;   The  proof  Is  by  contradiction  and  con- 
struction.  Assume  that  w^^  >  0  and  r^^  >  0  for  any  peri- 
od 1  gives  a  minimum  cost.   There  are  three  cases: 
Case  1   w.  -  r .  >  0 

Let  w'  =  W4  -  r.  and  r'  =  0 

1    ^    ^      1 
Then 


N  N 

Cj^w'-f  j^r'+(w'-r')*  J  m^  =  Cj^Cw^-rj^)-f^*  Of  (w^-ri) .  |  mj 

N 
=  Cj^Wj^-  c^rj^+Cwj^-r^^)  •  Z   m. 

N 
<  c^Wi-  firi+(wi-ri)-^Z^mj 


t-     J 
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since  c.>  f . . 

The  assumption  Is  contradicted  and  (1)  has 

been  constructed. 

Case  2   w.  -  r .  =  0 

Let  w'  =  0  and  r'  =  0 
1  1 

Then  certainly 


0  <  c^w.  -  f^r. 


since  Ci  >  fj^ 

The  assumption  Is  again  contradicted  and  (2) 
has  been  constructed. 
Case  3.  w.  -  rj^  <  0 

Let  w'  =  0  and  r*  =  r^  -  W4 

1  11^ 

Then 


N  N 

°i^'-^l^i"'^^l-^l)'j^i'^J  =  CjL-0-fi(ri-Wi)  +  (Wi-ri).^I^mj 


N 
^l*'l-^l^l"^^"i-'*l^'jSi'"J 

N 
T.   Ill  J 


<   c  j^Wj^-f  j^Tj^+Cwj^-r  j^)  •  .  ^  .ni 


since  c^    >  T^ 
The  contradiction  Is  complete  and  (3)  has 
been  constructed, 

QED 


:  -^e 
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Now  using  Standard  Transformation  1  as  In 
Section  II. 2,  constraint  sets  (P2.1)  and  (P2.2)  can  be 
combined  to  give  the  set  of  relations 

(w^-r^)  =  dj^  +  Yi  -  w^  -  v^ 

(w2-r2)  =  d2  -  d^  +  72  -  2yi  +  Vq 

(2.11) 

(w3-r3)  =  d^  -  d2  +  73  -2y2  +  yi 

(wi^-rjj)  -  dj^  -  dj^.i  +  Yn  -  ^yn~l    +  yN-2 

Note  here  that,  by  Lemma  2.9  and  w,  r  -  0,  If  the  right- 
hand  side  of  a  constraint  In  (2.11)  Is  positive,  then 
w.  >  0  and  r.  =  0.   If  It  Is  zero,  then  \i^   =  r^   =   0, 
If  It  Is  negative,  then  Wj^  =  0  and  r^  >  0. 

Solving  (2.11)  for  r  and  substituting  In  (P2) 
gives 


N 
P5:      mln  ^5=1  [piXi  +  h^Vi  +  f 'yi  +  (Ci-fj^)wj^]  +  M 


s.t . 

xi  ^  di+  y^  -  Vq 

(P5.1) 

x^  -  dj_  +  y^  -  yi_i  1  =  2,...,N 
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(P5.3) 


where 


x^  -  vi  =  d^  -  Vq 


(P5.2)     x^   +   Vi_i  -  v^  'd^  1  =  2,...,N-1 


^N  +  ^N-1  =  ^N  +  v^ 


^1  +  yi  -  yi-1  -  0  i  =  2,...,N 


^1  -  ^1  +  yi  -  (wq+Vq) 

(P5.4)     W2  -  d2  -  d-,^  +  y^  -  2yi  +  v^ 

w^  -  di  -  dj^_-L  +  y^  -  2yi_i  +  yi_2  i  =  3,..,N 
X,   1,    ^,    7   -    0        and  INTEGER 


M 


f^    =   f^    -    2fj^+;L    "^    ^1+2    +   n^i    -   "11+1        i    =   l,...,N-2 

^^^(^l-^i+l^^l    +    Vn   -^    Ji""!^!    -    ^l^o    ""    (f2-^l-"»l)^ 
Constraint    sets    (P5.1)   and    (P5.2)    and   the 


r  V 


•       *    / 
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first  two  sets  of  terms  In  the  objective  function  are 
Just  problem  P4  and  thus  are  equivalent  to  a  transporta- 
tion problem.   Writing  constraint  set  (P5.^)  separately 
with  the  last  set  of  terms  of  the  objective  function 
gives 


N 
P6:  mln  Zr   -      1    (c^-fi)  W4 

b    1=1   1   ^    1 


s .  t . 

w^  ^  di  +  y^  -  Wq 

w  ^  0  and  INTEGER 

which,  given  y.  Is  solvable  for  w  by  Inspection  since 
the  coefficients  In  the  objective  function  are  positive 
and  the  right-hand  sides  of  the  constraints  are  Integer. 
Furthermore,  It  Is  possible  that  the  right-hand  side  of 
some  constraint  1  Is  negative  or  zero.   Recalling  Lemma 
2.9,  It  follows  that  In  this  case  w^^  «=  0  and  r^^  equals 
the  absolute  value  of  the  right-hand  side.   So,  for 
each  possible  y  In  the  CIDPI  model,  It  Is  necessary  to 
consider  a  transportation  problem  and  a  trivial  linear 
program. 
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Essentially  the  same  branch  and  bound  enumer- 
ation scheme  can  be  followed  here  as  was  used  in  the 
capacity  increase  case.   Lemma  2.1  still  applies  as  a 
finite  bound  on  the  enumeration  and  (2.5)  is  still  a 
valid  set  of  global  upper  bounds. 

The  lower  bounds  on  y  are  of  a  slightly  dif- 
ferent form  than  (2.^)  and  can  be  gotten  from  (P5.3)  and 
y  -  0  as 

yi   -  v^  -  ^1 
(2.12)         yi  -  yi.i  -  d^     i  =  2,...,N 

II. 8  Properties  and  Algorithm  for  P5 

The  properties  of  P5  are  similar  to  those  of  P3. 
Because  of  the  new  lower  bounds  (2.12)  however,  most  of  the  re- 
sults must  be  restated  and  reproved. 

The  statement  and  proof  of  Lemma  2.2  hold  equiva- 
lent ly  for  P5  if  the  change  of  variables 

^i  =  (^1  "  ^1^      1=1,... ,N 
is  made.  To  arrive  at  sufficient  conditions  for  convexity  to 
hold  in  the  optimal  values  of  Z^,  the  remaining  results  must  be 
reiterated.  For  Lemmas  2.10,  2.11  and  2.12  and  Theorens  2.5  and 
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and  2.6,  when  given  a  set  of  y^,  1  =  1,...,J,  it  is  assumed 
that  the  remaining  y,_j_,    are  at  their  lower  bounds  for  k«l,...,N-J. 
Lemma  2.3  Is  restated  a^ 

Lemma    2.10      If   for   given  y^^,    1    «   1,...,J    the   lower 
bounds   of   yj+i^^  are   positive,    k  =   1,...,    N-J    then 
(wj+^-rj+i^)    =    0   for   k  =    2,...,    N-J . 

Proof:      For  any   ke{2, 3, . . . ,N-J}    from    (2.12) 

yj+k  -   yj+k-1   -  dj+k   >   0 
and 

yj+k-1   -   yj+k-2   -  ^J+k-l    >    0 
where   equality   holds  on  the   left.      Then   from    (2.11) 

^^j+k-^j+k)  =  (yj+k-yj+k+i+dj+k)-(yj+k-i-yj+k-2"^^j+k-i) 

«  0 

QED 

Next, Lemma  2.4  becomes 

Lemma  2.11   If  the  conditions  of  Lemma  2.10  hold  and 

y'  ■  yj"'"'^  ^°^  ne{l,2,3, , . .  >  then 


(wj-rj)'  =  (wj-rj)  +  n 


and  for  k  =  l,2,...,N-j  the  lower  bound  on  y'    is 

J+k 


,i  A  .      / 
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yj+k  +  " 

and 

and 

(wj+k-r-J+k)*    =0  k   =   2,...,N-J 

Proof:      From    (2.11) 

(wj-rj)'    =   dj    -   dj_-L    +    (yj+n)    -   2yj.i    +   yj.2 
=    Cwj-rj)    +   n 

For   k   =   1 

so  from   (2.11) 

^^j+l-^j+1^'    '^  ^J+1   ~  ^J    ■•"    (yj+l+i^)    -   2(yj+n)    +  yj_-L 
=    (wj+i-rj+i)    -  n 


Doing   induction  on   k,    for   k  =   2 

r»  ^    y' 

J+2  J+1 

J+2 
yj+2  +  n 


y:..  ^  y..,  -  dj+2 

=    (yj+i+n)   -  d 


I  , 


••      '. 
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Assume  the  relationship  holds  for  a  general  k. 
Then  for  k+1 


^1+k+l  "  ^.i+k  "  ^J+k+1 


J+k+1   -^J+k 

=  (yj  +  k-^n)  -  dj^j^^^ 


=  V      +  n 


Finally  by  Lemma  2.10 


QED 

The  following  is  an  additional  property  which  will  be 

needed. 

Lemma  2.12   If  the  conditions  of  Lemma  2.10  hold,  then 

^^'j+l-^J+l^'  -  0- 

Proof:   Prom  (2.11) 


since  y'    is  at  its  lower  bound  so  the  first  term  is 
zero,  and  the  second  term  must  be  non-negative  by  (2.12) 

QED 


Now  Theorem  2.1  can  be  restated  as 


I   ,, 
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Theorem  2.5   If  for  given  y^,  1  «  1,...,J,  the  lower 
bounds  on  y.^^^  are  positive  for  k  =  1,...,N-J,  then 
C(y,+n)  is  strictly  convex  In  n  for  nG{l,2,3, . . .} . 

Proof:   For  y»  =  Yj+n,  ne{l, 2,3, . . . }from  Lemmas 
2.11  and  2.12 

(wj-rj)'  »  (wj-rj)  +  n 
^^J+l-^J+1^*  "  (wj+1-^J+l)  -  n  ^  0 
(^^J+k-^j+k^*  "  °  k  -  2,...,N 

The  capacity  change  variables  in  period  J  could  start 
from  r^  >  0  and  go  to  w.  >  0  as  n  increases.   Assume 
there  is  an  n^^  such  that  r.  >  0  for  n  <  n-,  and  w^  >  0 
for  n  >  n,  ,   Then  the  change  in  capacity  cost  for  each 
increment  of  n  <  n-,  is 


'^1  ■  fj  *  "j  -  fj+i 


while  for  n  -  n-j^  the  cost  change  is 


i2  =  Oj  +  m_,  -  fj^^  >  0 


where  Ap  ^  0  follows  from  the  cost  assumption  of  P2 
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that 

^1  ^  J  A  >  ^J  +  k-j"^k     VJ  >  1.  1  -  1.....N-1 


Further  from  P2  it  was  assumed  that  c^  >  f^  giving  the 
convexity  result 


A-j_  <  A2 


QED 


Corollary  2.5.1   C(y  +n)  is  convex  in  n  for  ne{l,2, 

3  9  .  .  .  /  • 

Theorem  2.2  and  Corollary  2.2.1  still  hold 
since,  as  has  been  pointed  out,  the  production-inven- 
tory section  of  the  model  is  equivalent.   The  main 
result.  Theorem  2.3,  must  be  restated  as 
Theorem  2. 6   If  the  conditions  of  Theorem  2.5  hold, 
then  T(y^+n)  is  convex  in  n  for  ne{l , 2, 3 , . . . } . 

Proof:  The  proof  follows  from  Theorem  2.5  and 
Corollary  2.2.2  and  the  fact  that  the  sum  of  convex 
functions  is  convex. 

QED 
Corollary  2.3.1  is  unchanged.   The  above  properties  of 
the  P5  variables  and  convexity  in  the  optimal  values  of 
Zc  will,  with  a  few  exceptions,  permit  the  bounds  and 
optimality  tests  of  section  II. 4  to  be  used. 


f- 
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Algorithm  for  the  CIDPI  model 


The  bounds  and  optlmallty  tests  of  section 
11.^  hold  for  the  CIDPI  model  with  a  few  exceptions. 
Bound  3  Is  no  longer  valid  since  the  possibility  of 
decreasing  capacity  allows  a  w.  to  exceed  these  bounds 
and  not  violate  a.   For  Bound  h   the  lower  bounds  are 
now  found  from  (2.12).   Optlmallty  Test  3  can  no  longer 
use  Bound  3.   Finally,  since  Theorem  2.5  only  gives  the 
weaker  result  of  convexity  for  C(y),  Bound  8  has  to  be 
rewritten  to  ensure  that  the  test  Is  being  conducted 
In  the  area  where  the  slope  of  C(y)  is  positive. 
Bound  8_^  For  given  y^,  1  =  1,...,J  if  y^^  >  0,1  =  J+1, 
...,N,  C(yj)  -  C(yj+1),  and  C(yj+1)  ^  a   then  y^  Is 
fathomed  for  the  preceding  set  of  yj^. 

The  basic  algorithm  of  section  II. M  requires 
some  changes  for  the  CIDPI  model.   The  lower  bounds  In 
Steps  3  and  4b  are  now  found  from  (2.12)  and  the  Initial 
value  for  y,  in  Step  3  Is 

y-j_  =  max  [(vQ-d^),  0] 

The  first  part  of  Step  5  must  be  rewritten  to  reflect 
Bound  8'  as  follows:  With  PNT  =  J  determine  capacity 
cost  a  =  C(y.).   If  a  >  a  and  either  J  =  N  or  J  >  N, 


-',..-..,     J 
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y^    >   V    for  1  »  J+1,...,N  with  C(yj)  -  C(yj-l)  go  to 
Step  4a  with  PNT  »  J  (Bounds  8'  and  9) The  remain- 
der of  Step  5  Is  the  same.   Finally,  the  update  of 
Bound  3  in  Step  9  can  be  dropped. 

The  remarks  of  section  II. 5  pertaining  to 
alternate  optimal  solutions  are  relevant  here  also. 
Finally,  note  that  lower  bounds  (2.12)  will  serve  to 
increase  the  feasible  space  for  y  and  thus  longer 
average  running  times  should  result  for  the  CIDPI  model, 

II. 9  Cases  of  Monotone  Demands  with  Special  Cost 
Structures  for  the  CIDPI  Model 

As  with  the  CPI  model,  the  efficiency  of  the 

algorithm  for  the  CIDPI  model  is  increased  greatly  when 

monotonicity  exists  in  demand.   The  following  hold  for 
this  section 

a)  The  cost  coefficients  of  P2  (Chapter  I) 
are  monotone  non-increasing  with  time. 

b)  The  cost  assumptions  of  P2  hold. 

Monotone  non-decreasing  demand  case 

Lemma  2.13   If  assumptions  (a)  and  (b)  of  this  section 
hold,  demand  is  monotone  non-decreasing  and 
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<  N 
Nw  +  V   -   Id 
o    o    i=i 


i 


then  there  is  an  optimal  solution  with  yj^  =  0. 

Proof:   The  proof  follows  trivially  from  the  fact 
that  these  are  the  same  conditions  as  for  Theorem  2.4, 
Now,  In  addition  to  the  possibility  of  eliminating 
Idle  capacity  by  the  arguments  of  Theorem  2.4,  there  Is 
the  added  possibility  of  scrapping  It  and  recovering  a 
salvage  value,  so  the  results  hold  a  fortiori,   QED 

The  conclusions  of  section  11,6  which  follow 
from  Theorem  2.4  are  valid  here  for  the  result  above; 
l,e.  the  transportation  subproblem  becomes  trivial  and 
can  be  solved  by  Inspection  for  a  given  y,  and  a  much 
more  efficient  set  of  upper  bounds  can  be  Included 
implicitly  in  the  algorithm.   That  is.  Bound  lU  of  sec- 
tion II. 6  can  be  used  Instead  of  Bounds  6,  7,  8,  and  9. 
The  algorithm  changes  of  section  II. 6  remain  valid  with 
the  exception  that  the  new  lower  bounds  in  Step  4b'  are 
found  from  (2.12). 

Monotone  no n- Increasing  demand 

For  this  case  it  must  be  assumed  that 


^1  -  Vq  ^  d2 
in  order  that  the  monotonlclty  be  preserved. 
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Lemma  2 , m      If  demands  are  monotone  non-increasing, 
assumptions  (a)  and  (b)  hold  and  further 


k-1 
(c)        f J  +  I  m^   <    c^  Vk  >  J,  J  =  1,...,N-1 


then  w-]_  =  d-,  -  Wq  -  v  and  w.  =  0,  i  =  2,,..,N. 

Proof:   The  statement  for  w-j^  results  from  the  fact 
that  no  backordering  is  allowed  and  by  cost  assumptions 
(a)  and  (b)  and  monotone  demands,  it  can  never  pay  to 
build  extra  capacity  to  produce  early.   Assumption  (c) 
eliminates  the  possibility  of  scrapping  capacity  in  any 
period  to  the  extent  that  capacity  must  be  rebuilt 
later, which,  by  monotonicity  of  demand,  produces  the 
result  that  w^,  =  0,  i  =  2,...,N. 

QED 

Lemma  2.15   If  demands  are  monotone  non-increasing  and 
f  -  0  then  y^^  =  0. 

Proof:   The  proof  follows  trivially  from  the  fact 
that  demand  is  monotone  non-increasing  so,  with  f  -  0, 
any  idle  capacity  in  a  period  can  be  salvaged  for 
further  revenues, 

QED 

As  was  pointed  out  in  section  II. 6,  when 
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Yj^  =  0,  all  capacity  Is  used.   The  transportation  sub- 
problem  of  P5  thus  can  be  solved  by  inspection  for  a 
given  y,  speeding  up  the  CIDPI  algorithm  greatly.   The 
remaining  remarks  on  algorithm  changes  for  this  case 
are  the  same  as  for  the  preceding  case  of  monotone  non- 
decreasing  demands. 


Chapter  III 

In  this  chapter,  a  special  case  of  the  CIDPI 
model  will  be  presented  which  leads  to  a  formulation  of 
the  entire  model  as  a  network.   Other  extensions  of  the 
basic  model  are  then  presented  to  allow  piecewise- 
linear,  convex  production  costs  and  discontinuities  in 
capacity  acquisition. 

III.l  Network  formulation  of  the  CIDPI  model 

The  formulations  of  Chapter  II  required  a 
problem  decomposition  approach  to  provide  suitable  al- 
gorithms which  took  advantage  of  the  problem  structure. 
There  is  a  special  case  of  the  CIDPI  model,  however, 
where  the  formulation  can  be  simplified  to  that  of  a 
network  and  solved  with  a  transshipment  algorithm. 
This  occurs  when  the  cost  of  capacity  increase  is  equal 
to  the  salvage  value  in  any  period. 

Such  a  model  might  be  useful  when  there  is  a 
large  second-hand  market  in  the  units  of  capacity.   A 
firm  might  then  acquire  or  sell  units  with  equal  ease  in 
any  period  at  the  going  market  price.   Maintenance 
costs  then  reflect  the  investment  needed  to  keep  each 
unit  at  its  rated  capacity  for  possible  future  sale. 


J  ■■•    c 
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The  derivation  of  this  case  comes  from  recon- 
sideration of  the  CPI  model  PI.   If  w  -  0  is  dropped, 
allowing  W  to  be   unrestricted,  the  problem  would  still 
be  bounded  if  capacity  costs  Cj^,  i  =  1,...,N  were  mono- 
tone non-increasing,  since  from  constraint  (Pl.l)  and 
X  -  0  it  must  still  be  true  that 


1    > 

^ZqWj  ^0     i  =  1,...,N 


This  implies  that  each  w^  is  bounded  below  by  the  total 
capacity  built  up  to  period  i,  i.e. 


>  i-1 
w^  -  -  Z  w^     i  =  1, . . . ,N 

-L    j  =  o  ^ 


since  if  the  costs  are  monotone,  building  and  subsequent 
scrapping  of  capacity  can  never  result  in  a  profit. 
That  is,  with  c-j^  -  0  and  monotone  i  =  1,..,,N,  the 
right-hand  side  will  always  be  as  small  in  absolute 
value  as  possible.   In  fact,  as  a  limit,  it  could  never 
pay  to  build  more  capacity  in  any  period  than  what  is 
needed  to  produce  all  remaining  demand  in  that  period, 
since  the  excess  capacity  would  be  idle  and  can  never 
reduce  the  total  cost.   Thus,  with  the  assumption,  which 
will  hold  for  the  remainder  of  this  section,  that 

di  -  v„  i  0 
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a  set  of  lower  bounds  on  vJ  can  be  written  as 

wi  -  -Wo 

>  N 

>  N 

■^     i=J-l  ^ 

Now  with  ff  unrestricted  In  PI,  if  Standard 
Transformation  1  is  performed,  then  P3  no  longer  con- 
tains constraint  set  CP3.3).   The  integrality  of  y 
still  holds,  however.   To  transform  the  remainder  of  P3 
into  a  network  formulation,  slack  variables  s  -  0  must 
be  added  to  constraint  set  (P3.1).   Then  the  sum  of  set 
(P3«l)  is  subtracted  from  the  sum  of  set  (P3.2)  to  give 
a  redundant  constraint.   Finally,  multiplying  the  con- 
straints of  set  (P3.2)  by  -1,  the  transformed  problem 
can  be  written  as 


N 
P6:      min  Z^  «  Z  [pj^x^  +  hj^Vj^  +  c'y^]+  K 


s.  t. 

(P6.1)  xj^  -  y^  +  Si  -  d;L  -  Vq 

xi  +  y^^^   -  y^  +  s^  =  di     i  «  2,...,N 


/  ^ 
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-X,   +   Vt   =   -dn   +   V 

i      1        1      o 


(P6.2)      -x^  -  vi_i  +  v^  =  -di     i  =  2,...,N-1 

(P6.3)         Yn  -  J^si  =  0 

X,  y,  ^,  s  -  0  and   INTEGER 


where  the  constraint  (P6.3)  Is  redundant,  the  c'  and  K 
are  as  given  In  P3,  and  vj^  has  been  incorporated  ex- 
plicitly In  dj^. 

The  coefficient  matrix  associated  with  P6 
satisfies  the  sufficiency  conditions  for  total  unl- 
modularlty  as  described  by  Hu  (Hu  69).   This  guarantees 
that,  since  the  right-hand  sides  of  the  constraints  are 
Integer,  the  linear  program  associated  with  P6  without 
the  redundant  constraint  (P6.3)  will  terminate  Integer. 
P6  posses  more  useful  properties  however,  since  the 
problem  as  written  can  be  represented  by  the  network 
displayed  In  Figure  III-l.   Thus  the  problem  can  be 
solved  using  a  transshipment  algorithm  to  obtain  the 
minimum  cost  riow. 

The  only  evident  problems  with  the  network 
are  the  cycles  involving  y  and  s  which,  if  negative  in 
cost,  cause  unboundedness.   However,  since  flow  s^  is 
costless,  the  monotonlcity  assumption  on  the  cj  is 
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^1   O 


Figure  III-ll 
Network  for  model  P6 


The  figures  on  the  source  and  sink  arcs  represent 
the  arc  capacities.   For  the  remaining  arcs  the  figures 
represent  the  flow  variables  associated  with  the  arcs. 
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sufficient  to  prevent  this  as  seen  from  the  following 
lemma. 

Lemma  3« 1   If  capacity  costs  c'  are  monotone  non-in- 
creasing  for  i  »  1,...,N,  then 


N     > 


i^/;  -  0  .    J  =  1,...,N 


Proof:   It  is  obvious  that 


Cj^  +  mj^  -  0 


For  J  =  N-1 

(c|q_i-2cN+m|^.3_-ii^)  +  (cj^+mjq)  =  Ccj^_-j_-cjq)  +  ni^-l  *  ° 

and  by  induction  for  J  »  k,  1  -  k  <  N-1 

N  > 

i-k""!  '    ^''k-Ck+l)  +  mj,  -  0 

QED 

It  is  obvious  from  Figure  III-l  that  the 
only  cycles  are  of  the  form 

^N*  ^N-l*  •••»  yj »  Sj     ^    '   1,...,N 
which  have  non-negative  costs  by  Lemma  3.1.   This  con- 
firms the  fact  that,  since  P6  is  obtained  by  a  linear 
transformation  on  the  bounded  problem  PI,  it  too  must 
be  bounded.   Emphasis  must  also  be  given  to  the  fact 
that  the  capacity  change  costs  c',  incorporated  in  P6, 
are  completely  non-stationary  and  thus  admit  to  dis- 
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counting  or  any  other  aberrant  structure. 

Networks  and  the  general  CPI  and  CIDPI  models 

The  existence  of  unlraodularlty  In  the  special 
case  of  P6  permitted  use  of  a  network  algorithm  to  ob- 
tain an  optimal  Integer  solution.   The  general  models, 
PI  and  P2,  could  also  be  formulated  as  networks  but 
equivalence  of  the  networks  to  the  original  models  re- 
quires sets  of  constraints  external  to  the  networks. 
Moreover,  these  constraints  violate  the  unimodularity 
conditions,  as  shown  by  the  following  examples,  and 
thus  exclude  the  use  of  any  existing  network  algorithm 
to  gain  an  integer  solution.   These  examples  are  both 
linear  programs  having  non-integer  solutions  with  non- 
integer  objective  function  values  (eliminating  alterna- 
tive optima). 
Example  1^:   Example  for  the  CPI  model 

For  N»6,  3  -  (20,20,19,19,18,18),  p  «  (5,5,^,4,3, 
3).  h  »  (1,1,1,1,1,1),  IR  «  (3,3,3,3,2,2),  d  -  (3,7,9, 
13,18,21),  Wq  »  1,  Vo  »  1 

LP  optimum 

^1  =  2,  X2  =  7,  X3  =  xi|  =  X5  =  xg  =  15.25 

V3  =  6.25,  V4  =  8.5,  V5  =  5.75 
w-j^  =  1,  W2  =  5,  W3  =  8.25 
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Objective  Function  =  735.25 

Example  2^:   Example  for  the  CIDPI  model 

The  problem  is  as  in  Example  1  with  the  exception 
of  f  »  (8,8,7,7,6,6),  d  =  (3,5,8,14,9,7) 
LP  optimum 

XjL  =  X2  =  X3  »  xij  «  X5  =  7.6,  X5  »  7.0 

vi  =  5.6,  V2  =  8.2,  V3  =  7.8,  v^   =  1.4,  V5  =  0 

w-j^  =  6.6     r5  =  .6 

Objective  Function  «  457.2 

III. 2   Block  Capacity 

In  this  section,  CPI  and  CIDPI  models  will  be 
extended  to  handle  the  case  of  block  capacity.   This 
terminology  embraces  any  situation  where  capacity  can 
only  be  acquired  in  units  which  are  each  capable  of 
producing  large  numbers  of  units  of  demand  in  each 
period,  as  opposed  to  the  implicit  assumption  in  effect 
30  far,  that  each  unit  of  capacity  produced  a  single 
demand  unit  per  period.   Block  capacity,  then,  corres- 
ponds to  the  economists'  problem  of  "lumpiness"  in 
acquisition  which  typically  causes  serious  difficulty 
in  a  marginal  analysis  approach. 

In  the  remainder  of  this  section  it  is  assumed 
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that  In  any  period 

1  Capacity  Unit  »  k  Demand  Units 
where  k  is  a  positive  Integer  and  Is  stationary,  and 
that  Initial  capacity,  w  ,  Is  expressed  In  units  of 
capacity. 

CPI-'Block  capacity  model 

Working  first  on  the  CPI  model.  In  terms  of 
block  capacity,  (2.1)  can  be  rewritten  as 

kw^  -  Vl  -  d-L  -  C^Wq+Vq) 

k(Nw^+...+Wj^)  -  Vj^  -  dj^  +  ...  +  dj^  -  (Nkw^+v^) 

Following  the  same  development  as  In  Chapter  II,  and 
using  Standard  Transformation  1  to  solve  for  w  gives 


^  =  k  t^l  ■"  ^1  -  ^^V^o^^ 

W2  =  -  [y2  -  2y-L  +  d2  -  d;L  +  Vq] 

Ik 


Depending  on  the  source  of  the  demand,  management 
has  the  option  of  choosing  the  most  appropriate  unit  for 
demand,  e.g.  box,  pallet-load,  car-load,  etc.   The 
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and  the  expressions  for  w.  can  now  be  substituted  back 
into  PI  to  get  the  transformed  problem.   In  this  case, 
substitution  In  the  W  -  0  constraints  gives  somewhat 
expanded  results.   The  constraint  set  corresponding  to 
(P3.3)  now  looks  like 


-  [yj^  +  d^  -  (kw^+v^)]  L   0   and  INTEGER 


\    ^^U   '    ^^N-1  "^  ^N  -  ^N-1^  "  °    ^"^  INTEGER 

Taking  the  first  constraint  and  Ignoring  the  Integer 
requirement  gives 

^1  -  ^%  -^  ^o  -  ^1 
while  the  Integer  requirement  gives 


i  [y^  +  d^  -  (kw^+v^)]  =  n^ 


where   n.    is   a  non-negative  Integer.      But   this   Is   the 


same   as 


y,    =  n,k  -  d^    +   kw      +  v 
11  loo 


value  of  k  can  then  be  set  appropriately. 
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A  similar  development  for  period  J  gives 


and 


y  «  n^k  +  2y 


j-1  -  yj-2  -^  ^J-1  -  dj 


> 


for  n^  -  0  and  integer.   So,  in  general,  substitution 
in  w  -  0  gives  two  sets  of  constraints,  the  first  sim- 
ilar to  set  (P3.3) 


y  -  kw  +  V  -  d^ 
1     o    o    1 

(3.1)     y2  ~  2y^  +  ^1  -  ^2  ~^o 

yi  -  2y^_^  -  y^^2  +  ^i-l  -  ^i     i  •=  3,...,N 

y  -  0 


and  the  second,  a  set  of  integrality  constraints  on  y 


^1  =  ^1^  -  ^1  +  ^^o  ""  ^o 

y^  =  n_k  +  2y^  +  d,  -  do  -  V      n.  ^  0  and  INTEGER 
(3.2)  2    2      1    ^1    2    o      i 

i  -  1,...,N 
V  \^  ^  2y^.l  -'*yN.2  "   ^N-1  -  ^N 


The  effect  of  k,  then  is  to  superimpose  a  lattice  of 
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interval  k  on  the  space  of  y.   Up  until  now,  k  »  1  and 

so  the  enumeration  on  y  had  to  search  every  feasible 

integer  point  in  y  space.   As  k  becomes  larger,  however, 

the  number  of  enumerable  points  in  y  space  becomes 

smaller  for  a  given  d  and  the  problem  can  be  solved 

more  quickly. 

Doing  the  remaining  substitution  for  w  in  PI 

gives  the  following  form  of  the  CPI-Block  capacity 

model. 

P7:   min  Z_  =   I  [p.  x.+h.  v.  ]  +  -[  E  c'y  +K]  -  iiifllc-,w_ 


s.  t. 


(P3.1),  CP3.2),  (3.1),  and  (3.2) 


X,  V  -  0   and  INTEGER 


y  ^  0 


where  c*  and  K  are  as  in  P3.   The  basic  algorithm  is 
valid  for  P7  with  a  few  changes.   In  step  ^b  let 
yj  "  y.  +  k.   The  global  upper  bounds  are  changed 
slightly  by  the  parameter  k,  since  the  maximum  value 
for  Wj^  is  now 


1  N 
w,  =  -(  E  d^-v.)   -  w 
•"■    k  i«l  ^  ^      ^ 


I  ■     I-         1- 


J  , ,. 
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So  using  C3.2),  (2,5)  of  Chapter  II  becomes 

(3.3)     Yj  -  Jk«raax  w^  -  L  d^^  +  v^     J  »  1,...,N 


In  order  to  satisfy  (3.1) f  (3.2)  and  y  -  0,  the  lower 
bounds  on  y  are  first  obtained  from  (3.1).   If  y  -  0, 
then  (3.2)  is  also  satisfied.   If  y^  <  0,  1  -  i  ^  N, 
then  by  adding  ra.k  to  y .  for  successive  values  of 
n^  =  1,2,3,... >  until  y.  -  0,  a  feasible  lower  bound 
can  be  obtained  which  satisfies  (3.2), 


CIDPL-Block  capacity  model 

Now  for  the  CIDPI  model,  (2,11)  can  be  re- 
written in  terms  of  block  capacity  as 


1 
(wj^-r^^)  =  -[dj^  "•'  ^1  ~  ^^o"*"^©^-^ 

iv 

1 
(W2.r2)  -  -[d2  -  d-L  .  2y^   +  y^   +  v^] 


^^i-^i>  -  k^^i  -  ^i-1  +  yi-2  -  2yi-i  +  yi^ 


Substituting  for  r  -  0  in  P2  gives 


3,...,N 


^1  -  z^^i  +  yi  -  ^^^o-'^o^^  - 


and  INTEGER 


.    ..  .5 


^^d 
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^N ^^N  •  ^N-l  "*■  ^N"  ^^N-l  +  yN-2^  '   ^  ^^^   INTEGER 

Jv 


Following  a  development  similar  to  the  first  part  of 
this  section,  the  J^"  constraint  is  equivalent  to  the 
two  constraints 

^j  -  k^^j-^j-i^yj-2yj-i^yj-2) 

and 

yj  -  njk  +  dj_i  -  dj  +  2yj.-L  -  yj.2 


where  n^  is  any  integer,  but  now  n.  <  0  is  allowable 
since  the  right-hand  side  of  the  first  constraint  is 
no  longer  required  to  be  a  non-negative. 

In  general,  the  two  constraint  sets  which 
result  are 

^1  -  ^^^I'^yi-^  ^o-^^o)^ 

(3.M)  w^  -  j^C<i2-^i-2yi+y2+VQ] 

^i   '  ;:'^^l"^i-l'^yi~2yi-l'^yi-2^        ^   '    3,...,N 


and 


(3.5) 


^1  ""   "l^  -   d-,^  +    (IcWq+Vq) 

y2    «   n^)^   +    2y-j^   +    d^   -   d2   -   Vq 
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y      »   n^k  +    2y2    "   yi   +    ^2   ""    ^3  "l        INTEGER 

1    =    1,2,.. . ,N 

^N    =   '^N^  +    2yjj_i   -  yj^.2   +    %-l   "   "^u 

Substitution  for  r  in  the  remainder  of  P2  Just  repeats 
the  first  three  constraint  sets  of  P5  and  gives  the 
model 


N 
P8:      mln  Zg  =   E  [piXj^+hiVj^+f "y  +(Cj^-f ^)Wj^]  +  M' 


s.t. 


where 


(P5.1),  (P5.2),  C3.^),  and  (3.5) 
^,  7,  w  -  0   and  INTEGER 


f"  =  -f» 
i   k  1 


M'  =  -M  -  Aiizilf.w 
ir       k    1  o 


for  f»  and  M  from  P5. 
1 

Modifications  to  the  basic  algorithm  for  P8 
are  the  same  as  for  P7  except  that,  If  the  lower  bound 
for  any  y^  is  positive,  then  (n^^k)  should  be  subtracted 
from  y^  for  successive  values  n^  »  1,2,3,... »  until  the 
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smallest  y^^  -  0  is  obtained.   This  corresponds  to  the 

unrestricted  n  In  constraint  set  (3.5)  which  Is  caused 

by  the  possibility  of  having  r^^  >  0  In  any  period. 

It  should  be  noted  here  that  the  special 
cases  of  Chapter  II  are  valid  here  only  for  k  =  1. 
Another  point  worth  noting  Is  the  Increase  In  danger  of 
rounding  off  from  the  related  linear  programming  solu- 
tion as  k  Increases.   Not  only  are  the  unit  capacity 
costs  larger  for  large  k,  but  also  shifting  a  unit  from 
one  period  to  another  has  a  proportionally  greater 
effect  on  the  optimal  production-Inventory  schedule. 

Finally,  the  Increase  In  the  lattice  Interval 
with  k  leads  to  a  useful  scheme  when  deciding  on  a  capac- 
ity Increase  (and  decrease)  strategy.   If  several  sizes 
of  capacity  are  available,  then  the  algorithm  can  be  run 
for  each  value  of  k  with  the  running  time  probably  de- 
creasing for  each  larger  value  of  k.   This  allows  a  third 
dimension  of  optimization  for  management  with  the  cost 
of  gaining  the  added  Information  governed  by  the  algo- 
rithm running  time  for  the  smallest  value  of  k.   Note, 
however,  that  the  form  of  Standard  Transformation  1  does 
not  permit  k  to  be  non-stationary  within  a  model. 

III. 3   Convex  Production  Cost 

Up  until  now.  It  has  been  assumed  that  all 
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production  was  being  carried  out  by  one  shift  working 
regular  hours  with  unit  production  costs  p.   In  this 
section  the  CPI  and  CIDPI  models  with  block  capacity 
will  be  extended  to  allow  multi-shift  and  overtime  pro- 
duction.  Although  the  work  done  here  is  for  one  addi- 
tional  shift,  further  extensions  are  trivial,  as  will 
be  pointed  out  later. 

The  crucial  assumption  of  this  section  is  that 
the  unit  cost  of  production  increases  for  each  shift 
and  for  overtime.   In  the  manufacturing  context,  this 
assumption  is  not  restrictive  and  actually  models  real- 
ity, since  it  is  commonly  accepted  that  there  exist 
"shift  differentials"  in  manufacturing  costs  for  a 
second  and  third  shift.  These  can  be  ascribed,  in  part, 
to  the  increased  time  spent  on  maintenance  and  to 
greater  human  error  in  succeeding  shifts.   As  for  over- 
time production,  the  wage  differential  alone  is  suffi- 
cient to  Justify  the  assumption. 

Let 

X*-  second  shift  production  in  period  1 
p'-  second  shift  unit  production  cost  in 
period  1 
Production  capacity  has  been  defined  in  terms  of  regu- 
lar time  production.   It  is  assumed  that,  in  terms  of 
the  second  shift,  an  equal  amount  of  capacity  is 


.'   ') 


• .  r  • 
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available.   Thus 


<   1 

x'-   Iw.+w^     i»l,...,N 


Assume  that 


p^  <  p»    .  1  "  1,...,N 


Then  the  production  cost  function  for  period  1 


P^Cx^.x^)  =  p^Xj^  +  p^x^ 


Is  plecewlse-llnear  convex.   Similarly,  with  a  suitable 
rearrangement  of  terms,  total  production  cost 


P  =   L  P.(x.,x') 

1=1  ^   ^   1 


Is  plecewlse-llnear  convex. 

CPL-Blocjc  capacity-  Convex  production  cost  model 

First  the  extension  of  the  CPI  model  with 

N 
block  capacity  will  be  presented.   The  term  Z   p*x* 

1«1  1  1 

must  be  added  to  the  objective  function  of  P^^,  while 
the  following  constraints  are  added  to  the  constraint 
set: 
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1 

x'  -  k  L  w,  i  kw^     1  «=  1,...  ,N 
1     J»l  J      ° 

(PI. la) 

x'  -  0   and  INTEGER 


and  the  appropriate  block  capacity  coefficients  are  in- 
serted in  the  remainder  of  the  model.   Then  (Pl.l), 
(PI. la),  and  (PI. 2)  can  be  combined  to  give  a  set  of 
necessary  and  sufficient  conditions  for  feasibility 
which  take  into  account  two-shift  production. 


2kw^  -  Vt  -  di  -  v^  -  2kw_ 
1    J-    1    o      o 

2k(2w,+Wo)  -  v^  -  d,  +  d,  -  v^  -  Mkw 
1^2     12o      o 


>   N 
2k(Nw,  +  ..  .+w,,)  -  v^,  -   Z  d.  -  v^  -  2Nkw 
1      N     N    1=1  1    o       o 


Again  by  Standard  Transformation  1  a  solution  for  w  can 
be  obtained.   Then  substituting  into  w  -  0  of  PI  gives 
the  following  two  sets  of  constraints 


y.  -  v^  +  2kw  -  d. 
1    o      o    1 

(3.6)      y2  -  2y^  +  ^1  "  ^2  ~  ""o 

y^  -  2Vi^l   "   ^1-2  "^  ^i-1  ~  ^i     i  »  3,...,N 

and 


•v  t- 


r ,  ■ 
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y^  »  2n^k  -  d^  +  2kvi^   +  v^ 


72  -  2n2k  +  27-^  +  ^i  -  ^^2  "  ^o 


(3.7) 


> 


^  ^  ^     J.     ^     J       ^^^    INTEGER 

1  =  1,^,...,N 

^N  "  2^N^  "^  ^^N-l  -  yN-.2  "^  "^N-l  -  ^^N 


Prom  (3.7)  it  is  clear  that  the  lattice  Interval  on  the 
y  space  has  been  Increased  to  2k  with  the  Introduction 
of  a  second  shift. 

Completing  the  substitution  and  transformation 
from  PI  gives 


N 
P9:     min  Zq  =  I  [PiX^  +  p'x'  +  h^v^  +  cjy  ]  +  k' 


s.t . 


^1  -  i^^i  ■"  ^1  -  ''o^ 


(P9.1) 


X  -  -(d^  "^  ^1  ~  ^i-l^  ^  '  '^••••»N 


(P9.2) 


X.  <  i(d^  .  y^  .  v^) 
^I  -  i^^i  "   ^i  -  ^1-1^ 


i  -  2,...  ,N 


(P9.3)    x^  +  x^  +  v^^-^  -  ^i  '  ^i         ^  °  l.---.^ 

(3.6) 


t  ■•; 
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(3.7) 


X,  x',  V  -  0   and  INTEGER 


y>-   0 


where 


and 


c"  =  ir-c'     1  «  1,...,N 
1    2k  i  »    » 


for  c'  and  K  from  P3. 

1 

The  addition  of  the  extra  term  to  the  ob- 
jective function  and  the  extra  set  of  constraints 
modifies  the  transportation  problem  structure  of  the 
subproblem.   In  effect,  N  additional  sources  have  been 
added  and  the  number  of  arcs  has  been  doubled.   This 
development  might  make  it  beneficial  to  reformulate 
the  sub-problem  as  a  transshipment  problem  which  can 
be  solved  with  an  appropriate  minimum  cost  flow  algor- 
ithm (Ford-Fulkerson  62).  An  investigation  would  have 
to  be  made  of  the  tradeoff  between  the  faster  algorithm 
but  higher  network  growth  rate  of  the  transportation 
problem  against  the  slower  algorithm  but  lower  network 

growth  rate  of  the  transshipment  problem.  Other  algo- 
rithm changes  and  remarks  on  the  bounds  for  y  are 


;u        r 


.£   ■-.■     I 
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Identical  to  those  for  the  CPI  model  of  the  previous 
section. 

CIDPI-Block  capacity-  Convex  production  cost  model 

The  convex  production  cost  case  of  the  CIDPI 

model  with  block  capacity  can  be  obtained  easily  by  a 

similar  development.   In  P2 ,  the  term  l   p'x'  Is  added 

1=1  i  i 
to  the  objective  function,  the  constraints 


1         < 

(P2.1a)     x»  -  k  E  (w.-r.)  -  kw^     1  =  1,...,N 
1     1=1  J   J      o 


x'  ~  0    and  INTEGER 


are  added  to  the  constraint  set  and  the  appropriate 
block  capacity  coefficients  are  Inserted  In  the  remain- 
der of  the  model.   A  set  of  necessary  and  sufficient 
conditions  for  feasibility  Is  then  obtained  from  (P2.1), 
(P2.1a) ,  and  (P2.2)  as 

2k(w^-ri)  «  v^  ^  d^  -  Vq  -  2kw^ 

2k[2(w-L-r3_)  +  Cw2-r2)]  -  v^  -  d^  +  d2  -  v^  -  ^kw^ 


2k[NCw^-r^)  +  ..,  +  (Wj^-rj^)]  -  v^^  -  ^I  d^^  -  v^  -  2Nkw^ 


Using  Standard  Transformation  1  to  solve  for  r  and 


I    o  ■ 
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substituting  for  r  -  0  in  P2,  the  following  two  sets 
of  constraints,  analogous  to  (3.^)  and  (3.5) >  are  ob- 
tained. 

>   1   r 

w,  -  -k£d  +y^-v  -2kw  J 
1    211   o    o'' 


(3.8)     W2  -  ik[d2~d^-2y^+v^] 


w^  -  ;^t:d^-cii.i+y^-2yi_-L+yi.2]   ^  "  3,...,N 


and 


y  =  2n,k  -  d,  +  v  +  2kw 
1     1     1    o      o 

y  =  2n2k  +  2y-,  +  d-,-  dp  -  v 
(3.9)      2     ^      1    1   ^    o 

y^  -  2n3k  +  2y   -  y^  +  d^  -  d3    n^  INTEGER 

;  i  =  1.2,. ..,N 

y^  »  2n^k  +  2y^„i  -  y^^^   +  d^^^   -  d^ 


Again  note  that  the  lattice  interval  on  the  y  space 
has  been  increased  to  2k. 

Completing  the  substitution  for  r  in  P2,  the 
model  is  transformed  to  the  following: 

PIO:  min  Z.^   =  ^    [Pi x.+p' x'+h. v.+f" 'y  +( c.-f . )w. ]  +  M" 
■Lu   ^-2.     ■^^ii^'-'-ii   111 


V ,  ■^i-:    ■»-      \ 
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s.  t 


(PlO.l) 


(P10.2) 


^1   -   i^^l   ^   ^1  -   ^o^ 
X.    1    i(d^    .    y^    -    y^.,) 


1   =    2,... ,N 


1    «    2  , . . .  ,N 


(P10.3) 


x^  +   X'    +   Vi.-L  -  v^  =   d^ 


i  =    1,...,N 


CP10.4) 


^1  +  yi  -  ^o  -  0 


i  0 


1   =    2,...,N 


(3.8) 


(3.9) 


X ,    X »  ,   V ,   w  -   0        and   INTEGER 


y  i  0 


Where 


fit  t 

i 


2  k' 


i   =    1,... ,N 
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and 

where  f  and  M  are  from  P5. 

The  remarks  of  the  previous  section  relevant 
to  the  changes  In  the  algorithm  and  In  the  bounds  apply 
here  also.   As  noted  there,  the  subproblem  can  be  form- 
ulated as  either  a  transportation  problem  or  a  trans- 
shipment problem. 

As  has  been  mentioned,  the  techniques  of  this 
section  easily  extend  the  models  to  encompass  a  third 
shift  and  overtime  work.   The  capacity  available  for 
overtime  can  be  assumed  to  be  equal  to  that  for  one 
regular  time  shift.   This  presents  a  problem  since 
overtime  usually  Includes  both  time-and-one-half  and 
double-time  wage  periods.   Some  subjective  unit  pro- 
duction cost  would  have  to  be  chosen  to  represent  the 
entire  period  rather  than  differentiating  between  the 
two  cost  levels. 

Finally,  as  was  already  remarked,  an  add- 
ition of  another  shift  increased  the  lattice  interval 
on  the  y  space.   If  the  extension  is  continued,  then 
the  first  term  on  the  right-hand  sides  of  (3.6)  and 
(3.8)  becomes  3nk  with  a  third  shift,  and  4nk  with 
overtime,  increasing  the  interval  by  nk  each  time. 
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Chapter  IV 

The  basic  Idea  of  the  enumeration  algorithms 
presented  Is  similar  to  work  done  by  Murty  (Murty  68) 
and  Gray  (Gray  71).   A  discussion  of  these  related 
approaches  opens  this  chapter.   Also  Included  are  a 
short  discussion  of  the  restart  capability  of  the 
algorithms  and  a  presentation  of  computational  results. 
Finally,  some  conclusions  regarding  the  preceding  work 
and  suggestions  for  possible  extensions  are  given. 

IV. 1  The  Decomposition  Approach  to  Integer  Programming 

Since  the  algorithms  developed  In  Chapter  II 
use  the  Idea  of  decomposition.  It  would  be  useful  at 
this  point  to  put  them  In  perspective  relative  to  some 
similar  approaches  In  Integer  programming.   This  Is  not 
done  to  claim  extensions  to  existing  theory,  but  rather 
to  show  contrasts  and  similarities  occurlng  when  a 
common  approach  Is  used  on  different  problem  structures. 

One  approach  which  uses  the  properties  of 
duality  theory  to  develop  a  gradient  algorithm  for 
mixed  Integer  programming  was  given  by  Benders  CHu  69). 
A  drawback  to  the  algorithm  Is  that  efficient  use  of  the 
gradient  Information  requires  the  solution  of  a  series 
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of  Integer  programs  where  the  number  of  constraints  Is 
Increasing  by  one  at  each  Iteration,   The  algorithms 
of  Chapter  II  only  make  partial  use  of  the  available 
gradient  Information  In  the  convexity  properties  which 
are  developed.   Possible  further  use  of  this  Information 
was  not  Investigated  In  the  course  of  this  work. 

Other  approaches  more  closely  related  In 
principle  and  in  spirit  to  the  algorithms  of  this 
thesis  are  the  work  of  Murty  (Murty  68;  and  Gray  (Gray 
71).   Both  works  deal  with  an  exact  solution  to  the 
fixed  charge  problem  and,  as  pointed  out  by  Gray,  the 
two  have  mutually  enhancing  properties  and  could  pro- 
bably be  combined  to  give  a  powerful  algorithm. 

Murty  develops  a  decomposition  approach  which 
ranks  the  extreme  points  of  the  convex  hull  according 
to  the  value  of  the  variable  cost.   An  upper  bound  on 
variable  cost  is  given  by  the  use  of  a  minimum  value  for 
the  fixed  cost  in  a  manner  similar  to  that  which  gives 
an  upper  bound  on  capacity  cost  in  the  algorithms  of 
this  thesis.   The  variable  cost  upper  bound  can  then 
be  used  to  terminate  the  enumeration  of  extreme  points. 
The  algorithm's  efficiency  is  greatest  when  the  fixed 
costs  are  small  both  in  absolute  and  relative  magnitude 
with  respect  to  the  variable  costs  (Gray  71). 

Up  until  now,  no  mention  has  been  made  of  the 
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fixed  charge  aspects  of  PI  and  P2 ,   It  is  not  hard  to 

make  such  an  Interpretation  of  the  capacity  costs  how- 
ever, since  once  incurred,  they  are  spread  over  the 
units  produced  in  all  future  periods.   The  capacity 
cost  of  this  thesis  can  thus  be  thought  of  in  terms  of 
a  fixed  charge  although  the  problem  formulation  does 
not  make  this  explicit. 

In  light  of  the  above,  the  approach  of  this 
thesis  is  much  closer  to  that  of  Gray  than  of  Murty, 
since  Gray's  decomposition  approach  works  mainly  on 
tne  fixed  cost  section  of  the  problem  and  deals  with 
the  variable  cost  as  a  subproblem.   His  fixed  cost 
section  is  concerned  with  the  optimal  profile  of 
open  and  closed  routes  of  a  transportation  problem 
where  each  open  route  incurs  a  fixed  charge.   The  scheme 
used  is  enumeration  in  a  number  system  of  base  2" (where 
n  is  the  number  of  sinks)  which  takes  advantage  of  both 
the  transportation  problem  structure  and  the  0-1  prop- 
erty of  the  integer  variables.   In  Gray's  approach,  an 
upper  bound  is  developed  on  the  fixed  cost  which  is 
analogous  to  a  of  Chapter  II, 

The  major  difference  between  Gray's  (and 
Murty '3)  model  and  those  of  this  thesis  is  the  0-1 
nature  of  the  integer  variables.   Transformation  of  PI 
and  P2  to  0-1  variables  could  be  done,  but  at  the 
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expense  of  a  huge  Increase  In  the  number  of  variables. 

Gray  also  encounters  another  difficulty  not 
present  In  this  thesis.   He  has  no  means  of  ensuring 
that  each  node  of  the  enumeration  Is  feasible  to  the 
to  the  main  problem,  thus  expending  resources  In  enumer- 
ating, and  testing  for,  Infeaslble  nodes.   As  has  been 
mentioned  previously,  enumeration  on  y  assures  feasibil- 
ity at  every  node  In  the  algorithms  of  this  thesis. 

IV. 2  Restart  Capability  of  the  Algorithm 

The  algorithms  developed  have  a  feature  that 
enhances  their  flexibility  and  usefulness,  that  Is  the 
capability  of  starting  them  at  any  point  In  the  enumer- 
ation.  This  feature  Is  Inherent  In  the  method  of 
enumerating  on  y;  since  given  any  y'  and  the  value  of 
the  algorithm  pointer  (PNT)  plus  some  value  for  the 
current  bound  on  capacity  cost  (a;  and  the  current  op- 
timal cost,  the  algorithm  can  be  started  from  Step  5 
and  run  normally.   Even  If  a  and  current  optimal  cost 
are  set  to  Infinity,  the  next  execution  of  Steps  6 
through  9  will  establish  a  new  set  of  bounds,  though  If 
the  optimal  solution  precedes  the  point  at  which  the 
algorithm  Is  restarted,  the  new  bounds  may  not  be  as 
tight  as  those  which  would  have  resulted  from  the 
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optimal  solution.   If  the  optimal  solution  Is  given  by 
some  vector  after  y  it  will  be  found  whether  or  not 
Initial  bounds  are  provided. 

This  restart  capability,  when  combined  with 
the  linear  programming  bound,  provides  a  very  useful 
tool  to  management,  especially  when  a  large  problem 
which  could  be  expected  to  have  a  long  running  time  is 
encountered.   Instead  of  starting  the  algorithm  from  the 
lower  bounds  of  y,  a  y'  corresponding  to  a  w  slightly 
lower  than  any  rounded-off  solution  to  the  linear  pro- 
gram can  be  chosen  as  a  starting  point.   The  algorithm 
can  then  be  run  for  a  sufficient  period  of  time  to  in- 
spect all  w  which  could  have  been  rounded-off  from  the 
linear  programming  optimum.   An  optimal  solution  is 
assured  here  if  it  equals  Z*    .   If  none  such  is 
found,  then  a  management  decision  can  be  made  based  on 


the  difference  between  the  current  optimum  and  Z*  as 

LP 
to  whether  the  full  algorithm  should  be  run  in  an 

attempt  to  achieve  further  improvement. 

It  might  also  be  advantageous,  e.g.  because 

of  limited  computer  time,  to  run  the  algorithm  for 

relatively  short  periods  of  time.   If  no  terminal 

condition  has  been  reached  then,  again,  a  management 

decision  can  be  made  on  whether  to  attempt  to  improve 

the  current  optimal  solution  with  respect  to  Z*  . 

LP 
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If  so,  the  algorithm  can  be  started  right  where  it  last 
terminated  with  no  loss  of  efficiency. 

IV. 3   Employment  Smoothing 

The  relevance  of  the  literature  on  employment 
smoothing  was  mentioned  in  Chapter  I.   Now  that  the  de- 
velopment has  been  completed  in  the  context  of  capacity 
expansion,  a  reinterpretation  of  the  results  can  be  made 
to  handle  the  employment  smoothing  problem  lO'Malley 
et.al,  66,  Silver  67,  Lippman  et.al  67a  and  67b,  Lippman- 
Yuan  69)-.   Reference,  in  this  section,  will  be  to  the 
most  general  models  developed,  i.e.  those  of  section 
III. 3,  from  which  all  special  cases  can  follow  trivially. 

The  time  period  for  the  employment  smoothing 
problem  is  much  shorter  than  that  for  capacity  expan- 
sion by  the  simple  fact  that  hiring,  firing,  and  usually 
even  training,  labor  does  not  take  as  long  as  changing 
fixed  capacity.   It  is  assumed,  therefore,  that  fixed 
capital  capacity  is  never  a  constraint  on  the  size  of 
the  work  force. 

The  work  force  units  can  be  either  in  terms 
of  individuals,  or  groups  such  as  gangs,  crews,  etc. 
It  is  assumed  that  there  is  a  capacity  associated  with 
each  labor  unit  which  can  be  expressed  in  terras  of  the 
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consumptive  factor  bundles  d^.      This  capacity  could  be 
gotten  from  the  associated  fixed  capital  optimal  oper- 
ating level  capacity  or  by  a  separate  marginal  pro- 
ductivity level  argument  for  labor  as  was  done  for 
capital  In  Chapter  I.   It  Is  again  assumed  that  the 
production  function  exhibits  constant  returns  to 
scale  In  labor  and  the  consumptive  factor.   The  accom- 
panying assumption  Is  that  the  existing  labor  force  can 
be  used  at  any  level  up  to  Its  capacity  with  no  lost 
opportunity  cost  of  Idle  labor  other  than  wages. 

A  change  In  the  factors  of  capacity  will  have 
to  be  made.   Before,  the  capacity  definition  was  ex- 
pressed In  terms  of  one  shift  doing  regular  time  work. 
Now  capacity  Is  defined  In  terms  of  all  possible  regular 
time  output  which  can  be  achieved  by  the  existing  work 
force.   Thus  It  Is  assumed  that  the  shifts  are  filled 
sequentially.  I.e.  before  adding  second  shift  labor,  a 
full  first  shift  must  be  operating. 

The  CPI  model  P9  would  have  limited  appli- 
cation In  the  employment  smoothing  sense  because  It 
Ignores  capacity  reduction.   Its  simplifying  features 
relative  to  the  CIDPI  model  however,  make  Its  use  at- 
tractive when  the  circumstances  permit.   Two  Instances 
when  such  use  would  be  feasible  are  those  of  a  firm 
phasing  In  a  new  production  line  where  the  phase-In 
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period  is  relatively  long,  say  a  year;  and  that  of  a 
firm,  with  a  reduced  work  force  on  an  existing  line, 
facing  a  projected  growth  In  demand. 

Relnterpretatlon  of  P9  to  employment  smoothing 
leaves  the  Inventory  variables  v  and  their  coefficients 
h  unchanged.   The  production  variables  x  now  represent 
all  regular  time  production,  and  no  differentiation  In 
production  costs  p  among  shifts  Is  allowed.   Overtime 
production  is  x*  and  again  undifferentiated  overtime 
production  costs  p'  must  be  assigned.   The  capacity  ex- 
pansion variables  w  now  become  work  force  hiring  var- 
iables and  the  associated  costs  c  can  Include  all  the 
normal  costs  associated  with  hiring  personnel,  i.e.  ad- 
vertising, personnel  staff  costs,  training,  etc.   Wages 
now  take  the  place  of  maintenance  costs  for  m. 

As  the  reader  has  probably  noticed,  an  inter- 
pretational  difficulty  exists  concerning  the  capacity 
ascribed  to  overtime  production  in  P9.   The  assumption 
for  that  model  was  that  overtime  capacity  was  approx- 
imately equal  to  the  capacity  of  a  regular  time  shift. 
What  is  necessary  now,  if  Standard  Transformation  1  is 
to  be  used,  is  the  seemingly  ludicrous  assumption  that 
overtime  capacity  is  equal  to  all  regular  time  capacity. 
With  a  time  span  of  one  to  one  and  one-half  years,  how- 
ever, discounting  will  not  be  significant  so  that  over- 
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time  production  costs  will  be  much  higher  than  any  reg- 
ular time  costs  over  the  period.   Furthermore,  hiring 
costs  are  now  In  the  same  order  of  magnitude  as  pro- 
duction costs  so  that,  within  a  reasonable  range, 
hiring  additional  regular  time  labor  should  be  cheaper 
than  working  any  significant  amount  of  overtime.   A 
reasonable  approach  to  the  problem  would  be  to  first 
solve  It  for  regular  time  only,  and  then  try  a  solution 
permitting  overtime.   If  the  overtime  assumption  causes 
difficulty,  then  either  a  subjective  modification  of 
costs  can  be  made  and  the  problem  re-solved  or  the 
approach  has  to  be  abandoned.   For  the  majority  of 
cases  It  would  seem  that  the  approach  Is  useful. 

For  the  more  realistic  CIDPI  model  PIO,  the 
only  additional  Interpretation  necessary  Is  that  the  r 
now  represent  the  firing  variables  and  the  f  the  firing 
costs.   The  same  difficulties  exist  for  overtime  capaci- 
ty and  the  previous  remarks  apply  equlvalently  here. 

An  argument  can  be  made  for  the  equivalence 
of  hiring  and  firing  costs  In  any  period.   Certainly 
what  Is  Included  a».d  excluded  from  each  cost  Is  sub- 
jective.  Assuming  that  equality  exists  and,  further, 
that  the  c  are  monotone  non-lncreaslng  with  time,  then 
the  special  case  of  the  network  formulation  Is  appli- 
cable to  the  employment  smoothing  problem  with  the 
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preceding  interpretations  on  the  relevant  variables 
where  in  addition  w.  <  0  corresponds  to  the  firing 
variable  r.  >  0  of  PIO.   Moreover,  use  of  the  network 
excludes  the  possibility  of  overtime  and  requires  that 
k  =  1.   No  obvious  extension  of  the  network  model  to 
remove  these  restrictions  exists. 

Finally,  the  assignment  of  a  value  for  k  in 
P9  and  PIO  will  undoubtedly  be  more  difficult  in  terms 
of  labor  than  it  was  in  terms  of  fixed  capital.   On  the 
other  hand,  the  sensitivity  of  the  models  to  k  should 
not  be  significant  for  labor  because  the  hiring  and 
firing  costs  are  much  smaller  than  the  corresponding 
capacity  change  costs,  and  thus  the  impact  of  fluctu- 
ations in  the  labor  force  level  will  not  be  as  great. 

This  work,  then,  becomes  the  first  exact 
integer  approach  for  the  employment  smoothing  problem. 
The  algorithms  given  for  the  special  demand  cases  in 
sections  II. 6  and  II. 9  are  more  general  than  those  of 
Lippman  et.al.  (Lippman  et.al.  67b),  and  algorithms  to 
handle  general  demands  and  general  linear  costs  have 


-^As  was  mentioned  in  Chapter  I,  a  "round-off" 
approach  from  a  linear  programming  model  has  been  pre- 
sented (O'Malley  et.  al.  66). 
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been  given  in  addition.   Finally,  the  similarity  of 
the  two  areas  of  capacity  expansion  and  employment 
smoothing  certainly  suggests  that  a  combined  effort 
would  be  useful, 

IV. 4   Computational  Results 

The  algorithms  developed  in  this  thesis  were 
coded  in  Fortran  IV  and  run  on  the  IBM  709^  in  the 
Johns  Hopkins  University  computing  center.   The  test 
problems  used  were  set  up  by  the  author  using  non-random 
data.   Twelve  problems  were  generated  in  all,  each  one 
as  a  ten-period  problem.   When  run,  however,  each  prob- 
lem was  run  eight  times  using  the  first  k  periods  of 
data  each  time  for  k  =  :i,...,10.   Thus  e.g.  for  problem 
10,  the  three  period  problem  is  numbered  10.1, while  the 
ten  period  problem  is  numbered  10.8.   This  numbering 
system  is  followed  in  the  tables  of  data  which  follow. 

The  following  conditions  are  satisfied  by  all 
of  the  problems:  a)  Costs  are  monotone  non-increasing 
with  time.   In  addition,  the  costs  of  Problems  7,  ti,  9 
and  10  are  stationary,   b)   k  =  1  (1  Capacity  unit  =  1 
demand  unit),  c)  Only  one  shift  of  regular  time  work  is 
being  used. 

A  linear  program  was  solved  for  the  CPI  and 
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CIDPI  case  of  each  problem  (using  models  PI  and  P2). 
Because  of  the  "nearly"  unlmodular  structure  of  the 
general  models,  a  large  number  of  Integer  extreme  points 
probably  exist  in  each  convex  hull.   In  fact,  for  the 
166  linear  programs  solved  for  Problems  1  and  3  to  12, 
a  total  of  4  0  terminated  integer. 

It  had  been  thought  that  the  number  of  periods 
would  be  the  most  significant  factor  in  problem  running 
time  and  that  the  appropriate  curves  could  be  presented. 
Instead,  the  models  turn  out  to  be  sensitive  to  both  the 
cost  and  demand  structures  of  the  individual  problems, 
as  well  as  being  dependent  on  the  number  of  periods. 
The  computational  results  are,  therefore,  presented  in 
tabular  form  in  Tables  IV-1,  IV-2,  and  IV-3  at  the  end 
of  this  section. 

An  example  of  the  sensitivity  of  the  algorithms 
to  demand  and  cost  structure  is  given  by  Problems  9  and 
12  using  the  CPI  algorithm.   The  data  for  the  problems 
is  as  follows: 

Problem  9, 

d  =  (5,  7,  11,  6,  10,  14,  9,  13,  19,  16) 

c  =  (  30,  30,  ...,  30) 

ni  =  (  3.  3,  ...,  3) 

P  =  (  5,  5,  . . .,  5) 

H  =  (  1,  1,  ...,  1) 
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Problem  12 

d  -  (2,  3,  5,  4,  6,  8,  7,  9,  12,  10) 
c  -  (250,  2^0,  230,  220,  210,  200,  190, 

180,  170,  160) 
in  =  (3,  3,  .  . .,  3) 
P  =  (5,  5,  5.  4,  4,  ^,  4,  3,  3,  3) 
h  =  (1,  1,  ...,  1) 

For  t  =  running  time 

n  =  number  of  periods 
a  curve  fitting  program  was  run  for  a  polynomial 

t  =  an^ 

and  an  exponential 

bn 
t  =  ae 

fit  on  the  data.   For  Problem  9»  using  7  data  points  the 
results  were: 

Polynomial 

a  »  .00028     b  =  5.10 
Coefficient  of  correlation  =  .960 

Exponential 

a  =  .00635     b  =  .94 
Coefficient  of  correlation  =  .606 
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while  using  7  data  points  for  Problem  12  the  results 
were: 

Polynomial 

a  «  .00006     b  »  5.63 
Coefficient  of  correlation  =  .875 

Exponential 

a  »  .00573     b  =»  .87 
Coefficient  of  correlation  =  .973 

When  the  data  for  6  common  points  was  averaged  the 
following  resulted: 

Polynomial 

a  =  .00003     b  =  6.21 
Coefficient  of  correlation  =  . 56I 

Exponential 

a  =  .00339     b  =  .993 
Coefficient  of  correlation  =  .770 

Similar  results  were  obtained  using  the  data  from  the 
other  problems  for  both  the  CPI  and  CIDPI  cases. 
The  only  conclusion  which  can  be  drawn. 
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then.  Is  that  the  number  of  periods  involved  Is  a  sig- 
nificant factor  in  problem  running  time,  but  there 
exist  other  significant  factors.   Inspection  of  the 
tables  will  verify  the  claim  made  at  the  end  of  section 
II. 8  that  longer  running  times  should  result  for  the 
CIDPI  algorithm. 

Two  interesting  examples  remain  to  be  shown, 
the  first  demonstrating  the  acknowledged  fallacy  of 
rounding-off  the  linear  programming  solution,  and  the 
second  a  demonstration  of  the  power* of  the  restart  ca- 
pability discussed  in  section  IV. 2. 

Example  1:  For  the  CPI  case  of  Problem  12.2,  the  vari- 
able values  resulting  from  the  linear  and  integer  pro- 
gramming optima  respectively  are  (with  w^  =  1  and  Vq  =  1) 
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= 

0 

The  difference  between  building  all  capacity  in  period 
1  and  building  it  all  in  period  2  could  be  significant 
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to  management.   Furthermore,  the  production-Inventory 
policies  are  quite  different  in  each  case. 
Example  2:   The  CPI  algorithm  for  Problem  8,8  took 
27.067  seconds  to  terminate  where  termination  resulted 
from  the  linear  programming  bound.   The  optimal  capacity 
expansion  vector  for  the  linear  programming  solution 
was 

w  =  (1,  4,  3.875,  0,  ...,  0) 

The  CPI  algorithm  was  started  with  a  y  corresponding  to 
W-,  =  1,  W2  =  4,  w-^  =  3,  and  the  remaining  y^  through  y-,Q 
set  at  their  lower  bounds.   The  time  to  termination  by 
the  linear  programming  bound  from  this  starting  point 
was  5.9167  seconds. 

In  terms  of  core  usage  in  the  computer,  the 
algorithms  are  limited  only  by  the  requirements  of  the 
transportation  problem.   The  version  of  the  transporta- 
tion algorithm  used  requires  only  one  mxn  array  where 
for  the  algorithms  of  Chapter  II,  m  =  n  =  N,  while  if 
a  transportation  algorithm  is  used  for  the  convex  pro- 
duction cost  use  of  section  III. 3,  ni  =  2N,  and  n  =  N. 
The  enumeration  itself  is  performed  using  a  single 
vector  N+2  words  in  length.   The  growth  in  core  require- 
ments is  therefore  dominated  entirely  by  the  growth  of 
the  transportation  array. 
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Table  IV-1 

Computational  results  for  the  CPI  algorithm 
with  monotone  non-decreasing  demand 


Problem 

Periods 

(N) 

Optimal  Values  | 

Run  Tl 

.meCsec ) 

Step  8 
'Termination 

^J  umber 

^    LP 

IP 

LP 

IP 

1.2 

h 

385.0 

388 

.73 

.383 

No 

6.2 

4 

176.0 

179 

.77 

.20 

No 

7.2 

i^ 

211.0 

217 

.77 

.20 

No 

1.3 

5 

557.0 

558 

1.05 

1.533 

Ho 

6.3 

5 

250.0 

252 

1.03 

.40 

No 

7.3 

5 

309.0 

317 

1.05 

.317 

No 

1.^ 

6 

735.25 

736 

1.33 

.733 

Yes 

e,k 

6 

335.25 

336 

1.33 

.283 

Yes 

7.4 

6 

429.50 

437 

1.47 

.717 

No 

1.5 

7 

884.60 

885 

1.83 

1.383 

Yes 

6.5 

7 

418.40 

420 

1.82 

3.083 

No 

7.5 

7 

572.85 

578 

2.35 

2.467 

No 

1.6 

8 

1052.0 

1055» 

2.27 

8"^  min 

7.6 

8 

739.75 

741 

3.02 

11.383 

No 

6.7 

9 

583.60 

586 

3.50 

85.233 

No 

1.1 

9 

940.75 

942 

3.75 

75.350 

No 

6.8 

10 

668.25 

670* 

4.23 

11"^  min 

7.8 

10 

1173.0 

1174 

4.67 

474.50 

No 

*Thls  was  the  current  optimal  solution  when  the 
enumeration  was  stopped. 

^Step  8  uses   Z*   as  a  termination  rule  for  the 

enumeration. 
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Table  IV-2 


Computational  results  for  the  CPI  algorithm 
with  general  demand 


Problem 
Number 

Periods 
(N) 

Optimal  Values 

Run  Tlme(sec) 

Step  8 
Termination 

LP 

IP 

LP 

IP 

4.1 

3 

299.33 

310 

.40 

.283 

No 

9.1 

3 

373.0 

391 

.40 

.20 

No 

4.2 

4 

472.0 

478 

.73 

.883 

No 

8.2 

4 

271.0 

277 

.75 

.317 

No 

9.2 

4 

425.0 

445 

.72 

.383 

No 

12.2 

4 

668.75 

821 

.75 

.250 

No 

4.3 

5 

578.0 

581 

1.07 

2.533 

No 

8.3 

5 

377.0 

385 

1.07 

.70 

No 

9.3 

5 

513.0 

521 

1.10 

.517 

No 

10.3 

5 

232.0 

240 

1.07 

.433 

No 

12.3 

5 

855.40 

893 

1.07 

.450 

No 

4.4 

6 

636.0 

636 

1.43 

2.283 

Yes 

9.4 

6 

673.0 

678 

1.50 

1.650 

No 

10.4 

6 

313.50 

327 

1.52 

1.0 

No 

12.4 

6 

1099.0 

1185 

1.52 

1.050 

No 

8.5 

7 

568.14 

572 

2.35 

3.517 

No 

9.5 

7 

747.43 

750 

2.40 

2.10 

No 

10.5 

7 

346.0 

362 

2.30 

1.267 

No 

12.5 

7 

1245.86 

1268 

2.37 

1.850 

No 

Table    IV-2    (cont. ) 
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Number 

Periods 

LP 

IP 

LP 

IP 

Step  8 

^.6 

8 

781.0 

786 

2.87 

35.583 

No 

8.6 

8 

73^.0 

735 

3.05 

19.60 

No 

9.6 

8 

882.50 

884 

3.12 

3.867 

No 

10.6 

8 

M21.0 

422 

3.08 

1.817 

No 

12.6 

8 

1M49.0 

1534 

3.05 

4.233 

No 

8.7 

9 

843.25 

844 

3.83 

4.550 

Yes 

9.7 

9 

1101.25 

1102 

3.83 

16.433 

Yes 

10.7 

9 

523.0 

1 

525 

1 

3.80 

7.20 

No 

12.7 

9 

1738.89 

1844 

3.90 

21.867 

No 

M.8 

10 

921.75 

929» 

5.43 

lO"*"  min 

No 

8.8 

10 

1022.75 

1023 

5.62 

27.067 

Yes 

9.8 

10 

1261.50 

1263 

5.48 

287.550 

No 

10.8 

10 

603.25 

604 

5.45 

3.817 

Yes 

12.8 

10 

1917.50 

1965 

5.67 

41.50 

No 
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Table  IV- 3 

Computational  results  for  the  CIDPI  algorithm 
with  general  demand 


Problem 
Number 

Periods 
(N) 

Optimal  Values 

Run  Ti 

.meC  sec ) 

Step  b 
Termination 

LP 

IP 

LP 

IP 

4.1 

3 

299.33 

302 

.57 

1.067 

No 

9.1 

3 

373.0 

387 

.57 

.633 

No 

4.2 

4 

472.0 

472 

1.12 

10.80 

Yes 

8.2 

4 

271.0 

276 

1.13 

4.650 

No 

9.2 

4 

415.67 

429 

1.13 

9.217 

No 

12.2 

4 

668.75 

728 

1.15 

2.267 

No 

4.3 

5 

578.0 

58l» 

1.73 

9"^   min 

9.3 

5 

513.0 

521 

1.77 

255.350 

No 

10.3 

5 

232.0 

240 

1.73 

28.167 

No 

12.3 

5 

855.4 

868 

1.73 

38.350 

No 

9.4 

6 

673.0 

69  5» 

2.55 

8+  min 

10.4 

6 

313.5 

327 

2.53 

616.980 

No 

12.4 

6 

1099.0 

1113 

2.53 

880.0 

No 
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IV, 5   Conclusions 

The  algorithms  developed  In  this  thesis  take 
advantage  of  the  structure  Inherent  In  the  combined 
capacity  expansion  and  production-Inventory  problems 
when  no  backorderlng  Is  allowed.   The  decomposition 
approach  of  handling  the  fixed-charge  (capacity  ex- 
pansion) and  variable  cost  (product Ion- Inventory )  por- 
tions of  the  problem  separately  allowed  the  application 
of  special  techniques  for  each  part  which  could  then 
be  combined  Into  an  efficient  algorithm.   The  regular 
behavior  of  the  y  variables  In  the  fixed  cost  por- 
tion of  the  models  allowed  the  development  of  useful 
optimal  cost  convexity  properties  In  the  variable  cost 
section  and  In  the  overall  problem. 

The  main  problem  of  Integer  programming  al- 
gorithms has  not,  unfortunately,  been  avoided  In  this 
work.   That  Is,  these  algorithms  display  the  same 
unpredictability  In  running  time  displayed  by  all 
known  Integer  algorithms.   Furthermore,  computational 
experience  has  shown  that  the  algorithm  running  times 
are  sensitive  to  demand  and  cost  structures  as  well  as 
being  dependent  on  the  number  of  periods  Involved. 

The  flexibility  provided  by  the  restart 
capability  of  the  algorithms  is  one  of  their  most 
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attractive  features  from  a  management  point  of  view. 
The  feature  allows  for  efficient  use  of  the  linear 
programming  solution  Initially,  by  permitting  a 
search  among  all  rounded-off  values.   Available  com- 
puter time  can  then  be  used  efficiently  for  the  main 
algorithm  by  running  the  program  for  short  periods  of 
time,  since  no  Information  Is  lost  when  the  enumeration 
Is  stopped. 

The  economic  problem  of  lumplness  In  capacity 
acquisition  has  been  handled  rather  easily  in  the 
extensions  of  Chapter  III.   In  fact,  the  more  severe 
the  lumplness  is,  i.e.  the  higher  the  value  of  k,  the 
larger  the  intervals  in  y-space  where  the  enumeration 
is  performed. 

Finally,  the  observation  of  an  equivalence 
between  the  CPI  and  CIDPI  models  and  the  employment 
smoothing  literature  is  of  no  small  significance.   The 
amount  of  work  done  in  employment  smoothing  has  been 
fairly  extensive  and  unification  of  the  results  of 
these  two  areas  would  enhance  the  existing  results  in 
both. 

IV. 6   Extensions 

The  restrictive  assumptions  which  led  to  the 
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results  of  this  thesis  provide  for  some  obvious  ex- 
tensions.  Perhaps  the  one  with  the  greatest  usefulness 
In  the  management  context  would  be  a  provision  for 
backorderlng,  where  the  backorderlng  horizon  could  be 
restricted  to  M  periods,  for  M£{1 ,2 ,3 , . . . } .   Unfortu- 
nately, Incorporation  of  backorderlng  In  the  model  will 
destroy  the  structure  of  Standard  Transformation  1  and, 
In  addition.  If  a  decomposition  approach  Is  taken  to 
solve  the  problem,  the  size  of  the  enumeration  space 
would  be  greatly  enlarged.   On  the  surface  It  seems  as 
If  an  entirely  new  approach  must  be  developed  for  this 
extension. 

Another  useful  extension  would  be  the  develop- 
ment of  a  model  which  could  handle  concave  costs  of 
capacity  expansion.  I.e.  decreasing  marginal  cost  In 
the  addition  of  capacity  In  any  period.   The  optimal 
policy  results  of  Manne  and  Velnott  for  capacity  ex- 
pansion with  concave  costs  (Manne-Velnott  67)  might 
be  useful  here. 

As  was  mentioned  in  Chapter  I,  the  work  in 
capital  budgeting  is  Just  one  step  up  in  aggregation 
from  the  content  of  this  thesis,   A  simple  extension 
in  the  context  of  the  work  already  done  would  be  to 
have  a  capacity  expansion  budget  in  each  period.   This 
adds  N  budget  constraints  to  the  models  and  these 
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constraints  could  be  handled  easily  In  the  fixed  charge 
(capacity)  part  of  the  enumeration  algorithms. 

An  extension  to  many  production  lines  with  an 
overall  capacity    budget    in  each  period  would  not 
be  so  easy,  even  if  each  line  had  a  single  product. 
The  difficulty  will  spring  from  the  problem  of  opti- 
mally allocating  the  capacity  resources  among  lines  in 
a  period  as  well  as  optimizing  for  each  line  over  N 
periods.   A  very  large-scale  integer  program  will  cer- 
tainly result. 

A  smaller  scale,  but  equally  interesting, 
problem  would  be  an  extension  to  multiple  products  on 
the  same  production  line.   As  was  noted  in  Chapter  III, 
the  work  of  this  thesis  limits  k  to  a  single  value  for 
any  model  because  of  the  form  of  Standard  Transformation 
1,   Multiple  products  are  sure  to  imply  different  pro- 
duction rates  and,  hence,  different  values  of  k.   The 
challenge  appears  interesting. 


.      4.,    <  i    .. 
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